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FOREWORD 

Recent  requirements  for  increased  strength  and  service  life  of  machines  and 
structures  have  been  met  by  the  use  of  higher  strength  materials  and  new  fabrication 
and  joining  methods.  Simultaneously,  failures  due  to  fracture  have  increased  rela¬ 
tive  to  those  resulting  from  excessive  deformation.  Frequently  service  conditions 
are  such  that  low  temperature  brittle  fracture,  fatigue  fracture,  and  high  temper¬ 
ature  creep  rupture  must  be  considered  in  a  single  system.  National  concern  with 
increased  safety,  reliability,  and  cost  has  focused  attention  upon  these  problems. 

Methods  are  now  available  to  predict  both  fatigue  crack  initiation  life  and  crack 
propagation  life.  Paradoxically  the  materials  properties  required  for  long  fatigue 
crack  initiation  life  are  incompatible  with  the  requirements  of  high  fracture  toughness. 
Thus,  the  conflicting  design  approaches  and  requirements  placed  on  the  material  are 
confusing  and  often  impossible  to  satisfy. 

Numerous  publications  dealing  with  a  variety  of  fracture  problems  have  led  to 
many  new  and  useful  developments.  However,  the  synthesis  of  the  concepts  into 
methods  for  design,  testing  and  inspection  has  lagged. 

This  program  of  study  is  intended  to  contribute  to  the  integration,  correlation, 
and  organization  of  mechanics  and  materials  concepts  and  research  information  into 
a  form  that  will  permit  enlightened  decisions  to  be  made  regarding  fracture  control. 
Reports  are  in  preparation  in  three  categories: 

1.  Research  reports  designed  to  explore,  study  and  integrate  isolated 
and/or  conflicting  concepts  ard  methods  dealing  with  life  prediction, 

2.  Reports  to  introduce  and  summarize  the  state-of-the-art  concepts 
and  methods  in  particular  areas,  and 

3.  Example  problems  and  solutions  intended  to  illustrate  the  use  of 
these  concepts  in  decision  making. 

H.  T.  Corten 
Principal  Investigator 


ABSTRACT 


It  is  shown  in  this  study  that  linear  elastic  fracture  mechanics  techniques  do 
not  provide  a  consistent  geometry -independent  means  of  predicting  failure  in  micro- 
flawed  ductile  matrix  materials.  Several  aspects  of  microflawed  material  behavior 
during  fracture  testing  suggest  that  considerable  plastic  behavior  occurs  on  the  micro¬ 
scale  in  microflawed  materials.  It  is  suggested  that  limit  analysis  techniques  applied 
first  on  the  microscale  then  extended  to  the  macroscale  may  provide  a  consistent 
means  of  failure  prediction. 

Three-dimensional  limit  surfaces  are  developed  for  several  types  of  micro  - 
flawed  materials.  Comparison  of  these  limit  criteria  with  currently  available  the¬ 
ories  and  experimental  data  on  microflawed  materials  in  the  absence  of  macroscale 
notches  shows  that  the  proposed  theory  is  valid  for  smoo;h  material  failure  predictions. 

Subsequent  extension  of  the  microscale  limit  surfaces  to  continuum  macroscale 
structural  analysis  shows  that  the  proposed  theory  will  provide  an  upper  bound  failure 
prediction.  It  is  also  shown  that  linear  elastic  fracture  mechanics  methods  will  pro¬ 
vide  a  lower  bound  failure  prediction. 
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1.  INTRODUCTION 


In  the  past,  engineering  design  utilizing  microflawed  materials,  of  which  cast 
iron  is  a  prime  example,  has  been  based  on  elastic  strength  of  materials  concepts. 
This  amounts  to  neglecting  all  ductility  and  plastic  behavior  of  the  material  and 
usually  results  in  an  overdesigned  structure.  This  standard  design  procedure  for 
these  materials  has  come  about  because  of  the  different  properties  of  microflawed 
materials  in  tension  and  compression  and  because  of  an  inability  to  accurately  pre¬ 
dict  the  bending  behavior  of  these  materials. 

Typically,  microflawed  materials  have  different  moduli  of  elasticity  in  tension 
and  compression.  This  leads  to  a  noncentroidal  neutral  axis  and  non-linear  stress 
distributions  in  bending  situations  (1)*.  This  fact  manifests  itself  when  extreme 
fiber  failure  stresses  from  bending  tests,  calculated  using  elastic  formulas,  are 
compared  with  elastically  calculated  failure  stresses  from  tensile  tests.  Typically, 
the  elastically  calculated  bending  failure  stress  is  30  to  50  percent  higher  than  the 
tensile  failure  stress.  In  this  study  '‘failure"  is  defined  as  the  inability  to  support 
load  without  unlimited  deformation  or  fracture.  Failure  load  or  stress  is  defined 
as  the  maximum  load  or  stress  which  can  be  supported  by  a  part  without  unlimited 
deformation  or  fracture. 

The  increased  load  carrying  capacity  of  microflawed  materials  in  bending 
situations  is  usually  neglected  in  design.  Standard  procedure  is  to  simply  use  the 
tensile  yield  strength  as  the  limiting  extreme  fiber  stress  in  bending.  This  con¬ 
servative  approach  is  taken  in  design  because  of  an  inability  to  accurately  predict 
bending  failure  stresses. 

'■Numbers  in  parenthesis  refer  to  list  of  references. 
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The  problem  of  failure  load  or  failure  stress  prediction  in  microflawed  materials 
is  compounded  if  macroscale  holes  or  notches  are  present.  This  problem  has  not  been 
satisfactorily  approached  and  prediction  of  failure  stress  or  load  in  microflawed  ma  ¬ 
terials  under  various  loading  conditions  when  macroscale  notches  are  present  is  the 
subject  of  this  study.  For  purposes  of  the  following  discussion,  a  microflawed  ma¬ 
terial  is  defined  as  a  material  which  has  flaws,  either  holes  or  inclusions,  with  the 
smallest  dimension  not  smaller  than  10  microns.  This  definition  excludes  such  micro  ¬ 
scale  discontinuities  as  dislocations  but  includes  graphite  particles  in  cast  iron  and 
porosity  in  cast  metals.  Further,  discussion  will  be  limited  to  those  microflawed 
materials  which  have  a  ductile  matrix  capable  of  plastic  deformation.  This  rules 
out  materials  such  as  ceramics  in  which  plastic  behavior  is  negligible  or  non-existent. 

Over  the  past  10  to  ^5  years,  the  theoretical  tool  known  as  linear  elastic  fracture 
mechanics  (LEFM)  has  become  a  highly  developed  and  utilized  method  for  predicting 
failure  loads  in  materials  which  behave  in  a  brittle  fashion  when  a  notch  or  crack  is 
present.  The  plane  strain  fracture  toughness,  Kjc  is  the  major  parameter  used  to 
characterize  fracture  in  this  method.  is  thought  to  be  geometry  independent. 
Hence,  if  is  known,  LEFM  techniques  can,  in  principle,  be  used  to  predict 
failure  loads  of  a  notched  part  with  any  shape  macroscale  notch  under  any  loading 
condition.  This  tool  seems  to  be  ideally  suited  to  the  current  problem  and  is  the 
ob'.ious  method  to  apply.  Several  investigator  s  (2,3,4,  5,6)  have  attempted  to  use 
this  method  in  studying  failure  of  cast  irons  and  cast  steels.  The  combined  results 
of  these  studies  are  somewhat  inconsistent  as  will  be  discussed  in  Chapter  2.  In  an 
effort  to  clear  up  some  of  these  inconsistencies,  an  investigation  using  various  geom¬ 
etries  of  fracture  toughness  specimens  made  of  gray  cast  iron  was  undertaken.  The 
conclusion  of  this  LEFM  investigation  in  gray  cast  iron  was  disappointing  in  that 
accurate  failure  stress  in  bending  situations  could  not  be  predicted  from  information 
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obtained  in  tensile  tests  and  vice  versa.  At  this  point,  two  choices  were  available. 
Either  a  modification  of  existing  linear  elastic  fracture  mechanics  could  be  developed 
or  a  different  theoretical  tool  could  be  applied.  During  the  fracture  toughness  testing 
of  the  current  study,  considerable  plasticity  on  the  microscale  was  noticed.  Because 
of  this,  the  latter  tack  was  taken  and  the  theory  of  plastic  limit  analysis  was  applied 
to  the  problem. 

The  theory  of  Emit  analysis  was  first  developed  for  a  continuum  by  Drucker, 
Prager,  and  Greenberg  (7).  Subsequently,  various  investigators  applied  this  tech¬ 
nique  to  such  things  as  plates  with  holes  (8,  9, 10),  soil  failure  (11, 12),  reinforced 
concrete  beams  (13),  and  fibrous  composite  materials  (14, 15).  In  this  latter  work 
on  composite  materials  McLaughlin  and  Batterman  (14, 15)  first  developed  limit  cri¬ 
teria  for  small,  representative  elements  of  larger  structures.  Since  these  limit  con¬ 
ditions  were  set  up  for  representative  elements,  they  can,  under  certain  circumstances* 
be  used  in  continuum  structural  analysis.  The  conditions  under  which  the  representa¬ 
tive  element  limit  conditions  can  be  used  on  the  continuum  structure  depend  on  the  con¬ 
straints  imposed  on  the  structure  and  the  size  of  the  representative  element  when  com¬ 
pared  with  the  total  structure  (16). 

In  the  current  work,  this  concept  of  developing  limit  criteria  for  representative 
volume  elements  (  1VE  s)  on  the  microscale  and  then  extending  these  criterxa  to  mac¬ 
roscale  structures  will  be  employed.  Chapter  3  includes  a  development  of  the  various 
limit  conditions  for  two  types  of  RVE's  as  well  as  a  comparison  of  the  proposed  limit 
criteria  with  other  failure  criteria  and  with  data  available  in  the  literature  on  smooth 
specimen  tests.  In  Chapter  4,  these  RVE  limit  criteria  will  be  used  in  macroscale 
limit  analysis  of  bars  in  tension  and  beams  in  bending. 

The  theoretical  results  of  this  study  are  valid  for  any  microflawed  material 
which  has  a  microscale  structure  similar  to  the  RVE:s  to  be  studied  and  having  a 
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ductile  matrix.  However,  since  a  good  example  of  these  types  of  materials  are  the 
cast  irons,  the  comparisons  of  experimental  and  theoretical  results  discussed  herein 
will  be  for  cast  irons.  The  proposed  theory  adequately  predicts  the  failure  stress  for 
smooth  bar  situations  under  various  loading  conditions.  Further,  the  microscale  limit 
surface  can  be  extended  to  the  macroscale  in  order' to  provide  a  failure  prediction  tech 
nique  for  several  of  the  notched  specimens  used  in  the  current  study.  It  will  be  shown 
that  the  result  of  this  macroscale  extension  of  microscale  limit  criteria  to  cases  in¬ 
volving  macroscale  notches  is  to  provide  an  upper  bound  to  failure  load  prediction. 

On  the  other  hand,  it  will  be  shown  that  the  linear  elastic  fracture  mechanics  method 
provides  a  lower  bound  to  actual  failure  loads. 
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2.  LINEAR  ELASTIC  FRACTURE  MECHANICS  APPLIED 
TO  A  TYPICAL  MICRO  FLAWED  MATERIAL 

2.1  Introduction 

In  this  chapter, currently  available  literature, in  which  linear  elastic  fracture 
mechanics  techniques  have  been  utilized  in  analysis  of  typical  microflawed  material 
failure,  will  be  reviewed.  The  materials  studied  in  this  literature  have  been  gray  and 
nodular  cast  iron.  Also  to  be  presented  in  this  chapter  are  linear  elastic  fracture 
mechanics  data  from  a  gray  cast  iron  study  which  was  undertaken  for  this  current 
work.  It  will  be  shown  that  linear  elastic  fracture  mechanics  parameters  cannot  be 
used  in  predicting  failure  of  a  notched  microflawed  material. 

In  the  interest  of  brevity  and  clarity,  discussion  of  the  meaning  of  the  fracture 
mechanics  terminology  and  details  of  fracture  mechanics  testing  have  been  put  into 
Appendix  A  and  B  respectively,  rather  than  in  the  body  of  Chapter  2.  Hence,  in 

what  follows  it  will  be  assumed  that  the  reader  is  familiar  with  the  contents  of 
Appendix  A  and  B. 

—  MawedMwna'lf631  E'aStiC  FraCtUre  Mechailics  Applied  to  Typical  Micro- 

In  1969,  Clover  and  Pollard  (3)  presented  data  in  which  the  plane  strain  frac¬ 
ture  toughness,  of  a  gray  cast  iron  approximating  the  specifications  of  the  Amer¬ 
ican  Society  for  Testing  and  Materials  (ASTM)  for  Class  30  cast  iron  were  reported. 
The  reported  KJc  for  as-cast  material  with  normal  flake  size  and  100%  pearlite  matrix 
was  17.  5  ksi  \[Tn.  as  measured  in  a  four  point  bend  test.  In  the  work  from  which  the 
Glover  and  Pollard  publication  was  derived,  Glover  (2)  showed  that  the  "apparent  Kjc" 
did  not  vary  with  specimen  thickness  down  to  0.0984  inches.  The  Glover  (2)  and 
Glover  and  Pollard  (3)  "apparent  KJc"  is  taken  to  mean  the  conditional  Kfc  or  K 
obtained  from  experimental  load -displacement  P-v  records  (Appendix  B).  K  is 
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equal  to  Kj  only  if  the  P-v  record  meets  certain  criteria  on  linearity.  These  cri¬ 
teria  are  set  up  to  operationally  define  -  the  load  at  2%  crack  extension. 

Brandt  (4)  conducted  an  investigation  which  included  the  determination  of  plane 
strain  fracture  toughness  values  for  ductile  iron,  gray  iron  and  .25%  cast  steel. 

Table  1  lists  Brandt's  results  for  ductile  and  gray  cast  iron.  Brandt  fust  tested 
notched  three  point  bend  bars  at  -50°  F  which  were  5  inches  deep,  20.  5  inches  long 
and  2.5  inches  thick.  At  this  temperature,  Brandt  found  that  the  P-v  records  of  all 
three  materials  studied  were  non-linear  and  the  ASTM  criteria  for  valid  Kjc  testing 
could  not  be  met.  Several  tests  were  then  conducted  at  -320° F  where  more  linear 
P-v  records  were  observed  and  the  ASTM  criteria  were  valid.  Brandt  concluded  that 
the  materials  tested  had  a  good  deal  more  toughness  than  would  be  expected  from  strict 
application  of  the  ASTM  size  criteria  for  valid  testing  (see  Appendix  B), 

Lazaridis,  et  al.  (5,6)  conducted  an  investigation  on  several  types  of  ductile  cast 
iron.  For  purposes  of  comparison  of  fracture  resistance  of  the  various  types  of  iron, 
the  fracture  toughness  of  single  edge  notch  tension  specimens  (3/4  inch  thick)  was 
used.  One  of  the  ductile  irons  studied  is  approximately  equal  in  yield  strength,  tensile 
strength  and  elongation  to  the  ductile  iron  tested  by  Brandt  (4).  Table  2  lists  the  re¬ 
sults  of  Lazaridis  et  al.  for  this  ductile  iron  at  -40°  F.  The  average  value  is  higher 
than  Brandt's  Kq  listed  in  Table  1  for  three  point  bend  tests  at  -50° F.  However,  at 
-50° F,  Brandt  found  three  point  bend  tests  were  not  valid  whereas  Lazaridis  et  al, 
obtained  valid  tests  at  -40° F  on  single  edge  notch  specimens. 

Thus,  tests  at  approximately  the  same  temperature  but  for  different  geometries 
give  conflicting  results  in  terms  of  Kq  and  validity  criterion.  It  is  also  interesting 
to  note  that  if  the  thickness  for  valid  toughness  testing  (Appendix  B)  is  applied  to  these 
tests,  then  Brandt's  tests  on  2.  5  inch  thick  bars  should  be  valid  whereas  the  tests  of 
Lazaridis  et  al.  on  3/4  inch  thick  plates  should  not  have  been  valid.  Further,  in  both 
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cases,  the  secant  corrected  loads,  used  operationally  to  define  crack  extension 
and  hence  Kq,  are  considerably  below  the  final  fracture  load.  Thus,  if  Kj  can  be 
found,  LEFM  techniques  could  still  not  be  used  to  predict  failure  since  only  the 
load  for  crack  extension  would  be  defined. 

It  appears  from  the  sketchy  data  available,  that  may  not  be  geometry  inde_ 
pendent  in  cast  iron.  Further,  the  non-linear  behavior  of  cast  iron  may  mean  that 
can  be  used  to  predict  crack  advance  but  it  may  not  be  useful  in  predicting  fail¬ 
ure  loads  for  notched  parts.  In  order  to  investigate  these  points  further,  an  inves¬ 
tigation  using  fracture  toughness  specimens  of  different  geometries  was  undertaken. 

2.3  Current  Investigation  of  Linear  Elastic  Fracture  Mechanics  Techniques  Applied 
to  Gray  Cast  Iron  —  — — 

In  the  current  study,  several  geometries  of  fracture  toughness  specimens  were 

utilized.  The  material  used  was  ASTM  Class  60  pearlitic  matrix  gray  cast  iron  with 

type  A-l  graphite.  Figure  lisa  typical  engineering  stress-strain  record  for  the 

material  in  tension  and  compression  (17).  Table  3  lists  the  chemical  composition 

of  the  material.  The  fracture  toughness  specimens  were  taken  from  portions  of  a 

twelve  foot  diameter  cylinder  with  a  3.22  inch  wall  thickness.  Table  4  lists  the  var  • 

iation  in  proportional  limit  stress,  .2%  yield  stress  and  ultimate  stress  at  the  inner 

third,  middle  third  and  outer  third  of  the  3.22  inch  thickness. 

A  series  of  three  point  bend  (3PB)  fracture  toughness  specimens  were  tested 

at  various  temperatures  in  both  the  fatigut  precracked  and  not  fatigue  precracked 

condition.  The  specimen  dimensions  are  shown  in  Fig.  2.  The  curvature  shown  in 

the  figure  is  exaggerated  and  comes  about  due  to  the  six  foot  radius  of  the  cast  cylinder 

This  curvature  will  not  alter  the  test  results  since  it  is  small  in  the  actual  specimen. 

The  maximum  load  Pmax>  was  recorded  for  each  test  and  the  stress  intensity  based 

0n  Pmax’  called  Kmax  values  were  calculated.  Table  5  lists  the  results  of  the  test 
for  all  cases.  In  these  tests,  as  with  the  others  in  this  study,  the  stress  intensity 
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equations  were  obtained  from  Ref.  18.  No  appreciable  difference  was  noticed  in  the 
Kmax  values  between  the  machined  notch  and  the  machined  notch  plus  fatigue  crack 
cases  with  approximately  the  same  crack  length.  For  this  reason,  all  further  testing 
was  conducted  on  machined  notched  specimens  only.  The  three  point  bend  tests  also 
revealed  that  the  maximum  load  and  K  x  did  not  change  significantly  with  increased 
temperature  up  to  400°F,  Therefore,  all  further  testing  was  conducted  at  room  tem¬ 
perature. 

Several  tests  were  conducted  on  a  compact  (CT)  specimen  geometry  with  di¬ 
mensions  as  shown  in  Fig.  3A,  These  specimens  were  cut  from  segments  of  the 
cylinder  as  shown  in  Fig.  3B.  The  maximum  load  and  corresponding  values 

for  these  tests  are  shown  in  Table  6,  Three  additional  compact  specimens  were 
tested  and  the  load  versus  crack  mouth  displacement  records  were  taken.  These 
records  are  shown  in  Fig.  4  along  with  a  table  showing  the  secant  corrected  load, 

Pq,  the  maximum  load,  Pmax.  the  Kq  value  and  a  comment  on  whether  or  not  the 
ASTM  criteria  for  P-v  records  has  been  met. 

Circumferential  notched  tensile  (CNT)  specimens  and  double  edge  crack  (DEC) 
specimens  were  also  tested  in  the  current  study.  Figure  5  and  6  show  the  geome¬ 
tries  tested.  Table  7  lists  the  pertinent  specimen  dimensions,  fracture  load  and 
K  values  for  the  circumferential  notch  tension  test.  Table  8  lists  the  same  infor 

Hid  A 

mation  for  the  double  edge  crack  plate  specimens. 

In  terms  of  K  values,  the  circumferential  notch  tensile  data  from  three 

ulcl  X 

different  size  specimens  have  much  lower  toughness  than  the  three  point  bend  data. 

Table  9  summarizes  pertinent  test  data  for  the  four  specimen  geometries  tested. 

As  the  component  of  load  attributable  to  bending  decreases  from  the  three  point  bend 

to  the  compact  tension  to  the  double  edge  crack  plate  and  finally  to  the  circumferential 

notch  tension,  the  maximum  load  and  Kmgx  are  continually  decreasing.  Further,  it  is 

seen  from  Table  9  that  the  valid  secant  corrected  data  from  the  compact  specim.  ns 

agree  with  the  K  values  from  the  two  types  of  notched  tension  tests, 
max 

I 

. 
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Several  pertinent  observations  were  made  during  the  current  study.  First, 
catastrophic  failure  did  not  occur  suddenly  in  the  three  point  bend  or  compact  spec¬ 
imen  testing.  Instead,  the  cracks  extended  gradually  below  maximum  load  and  as 
the  load  was  increased,  the  cracks  grew.  There  was  no  "big  bang"  characteristic 
of  fracture  toughness  testing  except  in  the  direct  tension  tests.  Also,  microscopic 
observations  at  the  crack  tip  in  the  compact  testing  revealed  plastic  action  at  the  tips 
of  graphite  flakes.  This  was  evidenced  by  dimpling  of  the  matrix  areas  ahead  of  the 
flakes.  Another  investigator  has  noticed  this  (19)  and  Fig.  7  is  a  photograph  of  this 
effect.  Finally,  application  of  the  ASTM  criteria  on  required  thickness  of  fracture 
toughness  specimens  indicated  that  all  tests  should  have  been  valid.  Using  Kj  =  18.  8 
from  the  valid  CT  test  and  a  yield  stress  for  the  middle  of  the  casting  of  45  ksi  a  thick 
ness  of  0.435  inches  is  required  for  valid  tests.  Since  all  specimens  were  tested  in 
1  inch  thickness,  they  all  should  have  yielded  valid  KIc. 


Conclusions  of  Linear  Elastic  Fracture  Mechanics  Testing  and  Microflawed 
Materials 


Two  things  are  apparent  from  the  available  literature  and  the  current  testing 
program.  First,  all  information  indicates  that  standard  ASTM  procedures  do  not 
provide  satisfactory  results  when  used  in  fracture  toughness  testing  in  cast  irons. 
Secondly  and  more  importantly,  it  is  not  possible,  using  LEFM  techniques,  to  obtain 
geometry  independent  parameters  with  which  to  characterize  failure  in  cast  irons. 
This  can  be  demonstrated  by  a  simple  example.  In  the  current  study,  the  Kj  from 
the  compact  tests  is  18.8  ksi\/Tn.  This  value  can  be  used  as  the  critical  stress  in¬ 
tensity  in  an  analysis  of  a  three  point  bend  test.  As  an  example,  the  specimen  dimen¬ 
sions  of  3PB-1  listed  in  Table  5  and  the  stress  intensity  solution  for  this  geometry 
given  in  Ref.  18  will  be  used.  Solution  of  the  stress  intensity  equation  for  the  load  P 
gives  Pmax  -  7.  55  kips.  Comparison  of  this  calculated  maximum  load  from  the  tensile 
fracture  toughness  with  the  actual  fracture  load  of  the  specimen,  which  is  18.  8  kips 
from  Table  5,  shows  a  considerable  discrepancy. 
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One  reason  that  linear  elastic  fracture  mechanics  does  not  yield  a  good  predic¬ 
tion  of  failure  in  bending  is  due  to  the  different  elastic  moduli  in  tension  versus  com¬ 
pression  in  microflawed  materials.  The  formulas  of  linear  elastic  fracture  mechanics 
for  bending  situations  have  been  derived  assuming  a  material  with  equal  moduli  in  ten¬ 
sion  and  compression.  This  is  evidenced  by  the  use  of  the  elastic  bending  formula, 
a  =  Mc/I  in  the  various  bending  stress  intensity  derivations.  In  a  microflawe a  ma  ¬ 
terial  such  as  cast  iron,  the  modulus  of  elasticity  in  tension  and  compression  are  not 
equal  as  shown  in  Fig.  1.  This  means  the  neutral  axis  of  a  beam  is  not  at  the  centroid. 
Further,  Draffin  and  Collins  (20)  have  shown  that  the  neutral  axis  shin  s  during  loading. 
Thust  even  if  the  technology  of  linear  elastic  fracture  mechanics  could  be  modified  to 
account  for  the  off-center  neutral  axis,  some  method  of  accounting  for  shifting  of  the 
neutral  axis  would  have  to  be  developed  if  the  failure  load  is  to  be  predicted  accurately. 

This  shifting  of  the  neutral  axis  in  microflawed  materials  is  most  likely  why 
Pmax  and  Kmax  values  of  the  current  study  show  a  decrease  from  three  point  bend 
to  compact  to  direct  tension.  As  the  component  of  load  due  to  bending  decreases,  the 
amount  of  material  subject  to  compression  decreased,  the  neutral  axis  location  be¬ 
comes  less  important,  and  the  fracture  load  decreases  because  of  the  lower  strength 
of  the  material  in  tension  versus  compression. 

Two  alternatives  are  apparent  for  the  accurate  failure  prediction  of  microflawed 
materials  with  macroscale  notches.  First,  the  concepts  of  linear  elastic  fracture 
mechanics  can  be  modified  to  account  for  the  neutral  axis  shifting.  Some  of  the  cur¬ 
rently  available  non-linear  fracture  mechanics  techniques  such  as  the  J-integral  or 
R -curve  methods  could  be  useful  in  this  regard.  However,  these  techniques  have  been 
developed  with  the  aim  of  determining  fracture  toughness  on  specimens  smaller  than 
required  according  to  ASTM  criteria.  Indications  from  the  literature  and  current 
study  indicate  that  the  bend  specimens  tested  more  than  met  the  size  requirements 
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of  ASTM.  Further,  even  if  a  Jjc  and  hence  a  Kjc  could  be  determined  from  the  P-v 
record  of  a  bend  test,  failure  loads  could  still  not  be  predicted.  The  same  problem 
of  using  Kjc  to  predict  results,  discussed  earlier,  also  arises  here:  namely,  Kjc 
can  predict  crack  extension  buc  not  failure. 

The  other  alternative,  and  the  tack  to  be  taken  in  the  remainder  o*  tins  woik, 
is  to  investigate  the  apparent  microscale  plasticity  in  cast  iron  fUi,  es  and  attempt 
to  employ  plasticity  theory  on  the  microscale.  Two  important  factors  indicate  that 
considerable  plastic  behavior  is  exhibited  in  the  matrix  material  of  cast  irons.  The 
first  is  the  fact  that  considerable  non-linearity  is  noted  in  both  stress-strain  records 
of  tensile  tests  and  P-v  records  from  fracture  toughness  tests.  Tie  other  factor  is 
that  microscale  plastic  action  has  been  observed  at  the  tips  of  graphite  flakes  in  both 
the  current  investigation  and  others.  For  these  reasons  then,  it  appears  that  a  study 
of  plasticity  as  applied  to  microflawed  ductile  matrix  materials  may  lead  to  a  useful 
technique  for  failure  prediction. 
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3.  FAILURE  SURFACES  FOR  M1CROFLAWED  MATERIALS 
BY  LIMIT  ANALYSIS  TECHNIQUES 

3. 1  Introduction 

Because  of  microplastic  action  between  flaws  which  has  been  observed  in  micro- 
flawed  materials,  it  is  hypothesized  that  failure  of  such  materials  does  not  occur  until 
the  plastic  limit  load  of  small  elements  of  material  is  reached.  Therefore,  limit  anal¬ 
ysis  of  representative  volume  cements  will  be  performed,  and  the  resulting  limit 
loads,  expressed  as  a  surface  in  principal  stress  space,  will  be  in  effect  equivalent 
to  a  continuum  failure  surface  for  the  material. 

The  steps  to  be  followed  in  the  subsequent  analysis  are  illustrated  in  Fig.  8. 

The  first  step  is  to  develop  a  model  consistent  with  the  real  material.  This  model 
must  be  representative  of  the  real  material  and  macroscale  structures  must  be  made 
up  by  stacking  the  representative  elements  in  various  arrangements.  The  model  will 
be  called  the  representative  volume  element  or  RVE.  If  a  failure  or  limit  surface 
can  be  developed  for  the  RVE,  this  surface,  with  certain  restrictions,  will  be  a  valid 
limit  criterion  for  use  in  the  limit  analysis  of  any  notched  or  unnotched  macroscale 
structure.  It  is  the  purpose  of  this  chapter  to  propose  and  develop  a  failure  theory 
for  microflawed  materials  based  on  limit  analysis  applied  on  the  microscale.  Later, 
these  microscale  limit  conditions  will  be  applied  to  the  macroscale.  The  reasons  for 
using  this  approach  are  to  reflect  observed  microscale  plasticity  and  provide  a  theory 
which  would  involve  few  material  dependent  properties  and  would  be  based  on  micro¬ 
scale  features  such  as  flaw  size  and  spacing. 

The  theory  of  limit  analysis  was  first  developed  for  a  continuum  by  Drucker, 

Prager  and  Greenburg  (7).  The  two  major  theorems  of  limit  analysis  were  concisely 
stated  by  Drucker  (21)  as  follows: 

■ 
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Lower  Bound  Theorem  The  body  will  not  collapse,  or  will  be  just  at  the  point 
of  collapse,  if  an  equilibrium  distribution  of  stress  can  be  found  which  balances  the 
appaed  load  and  is  everywhere  below  yield  or  at  yield. 

Upper  Bound  Theorem  The  body  will  collapse  if  there  is  any  compatible  pattern 
of  plastic  deformation  for  which  the  rate  at  which  the  external  forces  do  work  equals 
or  exceeds  the  rate  of  internal  dissipation. 

The  major  points  of  these  theorems  are: 

1.  The  statically  admissible  stress  field  of  the  lower  bound  need  not  be 
the  actual  one  in  the  structure. 

2.  i  he  kinematically  admissible  velocity  field  of  the  upper  bound  need 
not  be  related  to  the  actual  velocity  field  in  the  structure. 

3.  There  is  no  required  relation  betwe  n  the  stress  field  of  the  lower 
bound  and  the  velocity  field  of  the  ~pper  bound. 

These  limit  theorems  have  been  used  by  many  investigators  to  obtain  failure 
load  prediction  in  such  things  as  plates  with  holes  (8,  9,  10),  soils  (11,  12),  and 
recently,  composite  materials  (14,  15).  It  is  proposed  here  that  some  of  the  tech¬ 
niques  used  in  these  analyses  be  used  in  obtaining  failure  load  predictions  for  micro- 
flawed  materials. 

McLaughlin  and  Batterman  (14)  and  McLaughlin  (15),  have  used  limit  analysis 
of  representative  microscale  elements  of  two  dimensional  structures  to  obtain  limit 
loads  for  macroscale  structures.  The  materials  these  investigators  studied  were 
fibrous  composites.  Majumdar  and  McLaughlin  (16)  have  shown  that  under  certain 
conditions,  limit  criteria  obtained  for  general  three  -dimensional  microscale  elements, 
called  representative  volume  elements  (RVE),  could  be  used  as  limit  criteria  for  con¬ 
tinuum  analysis  of  macroscale  structures. 


The  requirements  necessary  for  extension  of  microscale  limit  criteria  to 
macroscale  limit  analysis  as  set  up  by  Majumdar  and  McLaughlin  (16)  are  twofold. 
First,  the  microscale  element  analyzed  must  be  small  compared  to  the  larger  struc¬ 
ture.  Secondly,  the  constraints  and  type  of  loading  on  the  large  structure  must  be 
such  that  no  constraints  on  displacement  or  stress  distr  ibutions  are  placed  on  the 
RVE's  in  critical  areas  of  the  struct ure.  For  the  RVE’s  and  notched  structures  to 
be  considered  in  this  work,  both  conditions  are  met  and  the  micro  to  macroscale 
extensions  of  RVE  limit  surfaces  are  valid. 

rumit  conditions  for  two  types  of  RVE  will  now  be  developed;  one  type  is  a 
model  of  a  microflawed  material  containing  spherical  voids  while  the  other  is  a  model 
of  a  material  containing  multi -directional  slits.  This  latter  case  is  an  attempt  to 
model  cast  iron.  In  all  cases,  the  matrix  material  is  assumed  to  be  highly  ductile 
and  follow  the  von  Mises  criterion  with  equal  limit  strengchs,  a  in  tension  and 
compression.  All  RVE’s  are  three  dimensional  cubes. 

The  analyses  to  follow  are  for  periodic  arrays  of  the  RVE.  Hence,  strictly 
speaking,  the  materials  studied  should  be  anisotropic  in  normal  stress  space.  How¬ 
ever,  it  will  be  assumed  that  the  equations  developed  are  valid  for  statistically  iso¬ 
tropic  materials  with  a  random  distribution  of  flaws  under  principal  stress.  There¬ 
fore,  throughout  the  remainder  of  this  work  normal  and  principal  stress  on  the  RVE 
will  be  used  interchangeably.  Further,  stress  states  which  are  other  than  principal 
stresses  can  be  obtained  from  the  equations  developed  here  by  using  the  standard 
transformation  procedures. 

One  important  aspect  of  the  application  of  limit  analysis  techniques  in  studying 
the  problem  at  hand  must  be  pointed  out  at  the  outset.  This  is  that  limit  analysis  was 
originally  derived  assuming  an  elastic -perfectly  plastic  material  behavior  with  unlim¬ 
ited  ductility.  McLaughlin  (22)  has  recently  shown  that  the  techniques  of  limit  analysis 
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are  still  valid  in  work  hardening  situations  but  the  requirements  of  unlimited  ductility 
still  remains.  In  some  of  the  cases  to  be  studied  here,  this  requirement  may  not  per¬ 
mit  the  strict  application  of  limit  analysis.  In  cases  where  unlimited  ductility  does 
not  exist,  McLaughlin  (22)  has  shown  that  the  application  of  limit  analysis  still  pro¬ 
vides  an  upper  bound  to  the  failure  load.  These  facts  will  be  utilized  in  Chapter  4. 


3.2  Limit  Analysis  Applied  to  RYE  of  Structure  with  a  Periodic  Array  of  Holes 

Figure  9  is  a  sketch  of  the  RVE  model  for  the  case  of  a  periodic  or  macro  - 
scopically  homogeneous  array  of  holes.  The  RVE  is  simply  a  cube  with  side  dimen¬ 
sions  of  unity  and  1/8  spheres  at  all  eight  corners.  The  spheres  are  merely  hoios; 
hence,  they  have  no  yield  surface  and  zero  strength  in  tension  and  compression. 

This  model  would  be  useful  for  a  highly  porous  metal  such  as  bad  casting  or  certain 
portions  of  a  sintered  powder  product. 

The  lower  bound  analysis  of  the  hole  case  is  straightforward.  Assume  that 
structural  members  extend  through  the  structure  from  one  face  to  its  opposite  face 
in  all  three  directions  as  shown  in  Fig.  10.  For  the  assumed  normal  stress  field, 
there  are  seven  limit  conditions.  In  those  parts  of  the  structure  where  uniaxial 
stresses  only  occur,  these  conditions  are 
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=  CF 
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my 


(1) 

(2) 

(3) 


where  the  subscripts  1,  2,  and  3  indicate  the  principal  stresses  m  the  x,  y,  and 
z  directions  of  the  cube.  For  those  parts  of  the  assumed  structure  where  biaxial 
stresses  occur,  the  limit  conditions  are 


-  CT, 
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Finally,  in  those  areas  where  all  three  stresses  occur,  the  limit  condition  is 


[<al  -  a2 >2  +  <a2  '  <t3)2  +  <03  •  al >2]} 


1/2  _ 


The  limit  conditions  in  Eqs.  1-7  are  in  terms  of  actual  stress  on  the  microscale. 

In  order  to  get  the  equations  in  terms  of  average  macroscale  RVE  stresses,  it  is 
noted  that 

P1  P2  P3 

a  =  j—  i  ^2  ~  a-'  CT3  =  A — 

Am  m  m 

where  P2>  and  P3  are  the  loads  in  the  1,  2,  and  3  directions.  Am  is  the  area  of 
an  RVE  face,  which  is  unity,  minus  the  projected  area  of  all  flaws  within  or  on  the  RVE. 
The  average  stresses  are  defined  as 


al  =  K 


‘2  Ar 


°3  Ar 


where  the  bar  indicates  average  over  the  RVE  face  and  is  the  total  surface  area 

of  one  RVE  face.  The  P^  and  P2  values  are  equal  in  both  cases  so  that 


o.  A„  .  =  cr.  A 
1  Tot  1  m 


—  m 

"  "  ATot  "i 


A^  is  simply  1  since  a  unit  cube  has  been  considered  for  the  model. 

Tot 

The  seven  limit  conditions  for  the  lower  bound  can  now  be  rewritten  in  terms  of 
average  stresses  on  the  RVE  as 


*****  *■* 
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cr  =  A  cr 
1  m  my 


cr0  =  A  cr 
2  m  my 
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(a22  -  S2  53  +  a2)  =  (A^  amy)2 

‘  ff3  51  +  ffl2)  =  <Am  amf 


i  [((7j  -  a/  +  <a2  -  a/  +  (<r3  -  <Jj)2]  -  (Am  ®my) 


(12a) 


(12b) 


(12c) 


In  loading  the  RVE  from  zero  load  under  any  combination  of  o  ,  and  stress, 
whichever  limit  condition  is  reached  first  will  govern  the  lower  limit  load  failure. 
Thus,  the  lower  bound  three  dimensional  limit  surface  is  shown  in  Fig.  11(a).  Under 
a  plane  stress  condition,  with  =  0,  the  plane  stress  limit  condition  for  the  lower 
bound  appears  as  shown  in  Fig.  11(b). 

In  order  to  obtain  an  upper  bound,  a  velocity  field  must  be  assumed  for  the 
deformation  of  the  RVE,  then  the  external  work  is  set  equal  to  the  internal  energy 
dissipated  in  performing  the  assumed  deformation.  For  the  case  at  hand,  a  homo¬ 
geneous  deformation  of  the  matrix  will  be  assumed.  Hence,  the  internal  energy  dis¬ 
sipation  is  given  by 


e  .  dV 

ij  m 


(13) 


where  the  subscripts  i  and  j  are  the  standard  index  notation  with  a  range  of  three 

and  double  indices  indicate  summation  according  to  the  standard  summation  convention. 

V  is  the  volume  of  the  matrix  material  only.  It  is  noted  that  since  the  total  RVE 
m 
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volume  is  unity,  V  is  the  volume  fraction  of  matrix  or  one  minus  the  volume  frac- 
J  m 


tion  of  voids.  Since  the  model  axes  are  principal  axes  and  the  displacement  rates 
are  linear  functions  of  the  strain  rates,  Eq.  13  becomes 


Vm  <al  'l  +  a2  <2  +  ct3  V 


(14) 


where  cr^,  a ^  and  are  the  stresses  at  the  point  on  the  matrix  limit  surface  where 
the  normal  vector  has  components  proportional  to  an<^  Cg  in  the  o^,  a  and 

<t  directions,  respectively.  Letting  the  displacement  rates  in  each  of  the  three  prin¬ 
cipal  directions  of  the  RYE  be  u,  v  and  w,  the  external  work  done  is 
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where  the  displacement  rates  are  linear  functions  of  the  strain  rates,  that  is 

ii  -  X  e. 


1 


v  =  Y  e2 

w  =  z 


(16) 


Using  Eq.  16  in  Eq.  15  and  evaluating  the  limits,  the  external  work  becomes 


al  el  +  °2  c2  +  a3  e3 


(17) 


Setting  the  external  work  equal  to  the  internal  energy  gives 

5i  h  +  52  H  +  “3  '3  '  Vm  <al  +  a2i2  +  a3  k3> 


(18) 


Since  the  e  terms  are  completely  arbitrary,  the  coefficients  of  the  e  terms  can  be 
set  equal  giving 


.'•tJUj-U.Jj,. 


19 


cr  a  and  cj„  are  combinations  of  stress  which  satisfy  the  von  Mises  limit  condition. 

1  Z  O 

Hence,  in  terms  of  average  stresses,  the  upper  bound  limit  criteria  for  the  RVE  under 
the  assumed  deformation  becomes 


i  [OTj  -  a2)2  +  (a2  -  S3)2  +(a3  -  a/]172  -  (Vm  %y): 


In  a  ,  and  o,  space,  this  is  simply  the  von  Mises  cylinder. 

1  z  o 

It  will  be  noted  that  the  upper  and  lower  bounds  are  considerably  far  apart  in 
cases  of  triaxial  tension  or  compression.  The  lower  bound  closes  off  in  these  octants 
and  predicts  failure  under  hydrostatic  conditions  whereas  the  upper  bound  is  not  closed 
off  in  these  octants.  In  triaxial  stress  situations  involving  one  principal  stress  which 
is  of  opposite  sign  to  the  other  two,  the  lower  and  upper  bound  are  concentric  von  Mises 
cylinders.  The  difference  between  these  cylinders  depends  on  the  ratio  of  Vm  to  A^. 

In  order  to  close  off  the  upper  bound  in  triaxial  tension  or  compression,  a  ve¬ 
locity  field  must  be  developed  for  the  RVE  which  would  allow  deformation  of  opposite 
faces  without  requiring  displacement  of  faces  perpendicular  to  the  ones  displaced. 

This  is  necessary  in  order  to  obtain  limit  conditions  giving  flat  plane  cut  off  surfaces 
to  the  von  Mises  upper  bound  cylinder.  This  exercise  was  not  carried  out  in  this  work, 
and  further  consideration  of  this  RVE  will  be  restricted  to  plane  stress  conditions. 

In  plane  stress,  Eq.  20  describes  a  von  Mises  ellipse.  Comparison  of  the 

lower  and  upper  bound  surfaces  in  plane  stress  (o^  =  0),  Fig.  11(b)  conditions 

shows  that  the  ratio  of  lower  to  upper  bounds  in  biaxial  tension-compression  regions 

is  V  divided  by  A  .  Since  the  assumed  model  is  for  spherical  holes, 
m  :  m 
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Vm  1  -  4/3  7r  r3 

A  t  2 

m  1  “  7r  r 

where  r  is  the  radius  of  the  spherical  hole.  Hence,  for  small  r,  the  ^m/Am  ratio 
is  close  to  unity  and  the  bounds  are  close.  Sphere  radii  between  0.0  and  0.25  give 
upper  and  lower  bounds  that  are  within  13%.  The  maximum  radius  that  a  sphere  in 
the  unit  cell  can  have  is  0.  50. 

For  actual  materials  approximating  the  model  used,  upper  and  lower  bounds 
can  be  determined  by  estimating  the  flaw  size  and  spacing  using  a  photomicrograph 
of  the  structure  and  hence  determine  the  radius  of  the  flaws  to  use  in  the  RVE,  the 
load  carrying  area,  and  the  volume  of  matrix.  As  an  example  of  this  technique, 
assume  that  the  structure  of  a  hypothetical  microflawed  material  is  made  up  of  ran¬ 
domly  spaced,  uniform  size  holes  as  shown  in  Fig.  12(a)  (23).  If  a  plane  AB  is 
passed  randomly  through  this  structure,  the  plane  would  appear  as  in  Fig.  12(b). 

The  principal  values  which  must  be  determined  to  utilize  the  proposed  analysis  are 
the  flaw  radius  and  spacing.  The  chances  that  a  given  plane  passed  through  the  micro- 
structure  will  cut  any  flaw  at  the  equator  are  slim;  nevertheless,  the  largest  radius 
observed  in  the  plane  under  study  can  be  reasonably  used  as  the  flaw  radius  for  the 
model.  The  center-to-center  spacing  of  flaws  in  the  thre-  dimensional  structure  can 
be  determined  from  measurements  on  the  random  plane  (23).  This  is  accomplished 
by  drawing  random  test  lines  on  the  plane  and  determining  the  number  of  particles 
intercepted  per  unit  length  of  these  test  lines.  The  reciprocal  of  this  value  is  the 
three-dimensional  center-to-center  spacing.  For  the  purpose  of  illustration,  suppose 
that  in  the  hypothetical  structure,  the  largest  flaw  on  a  plane  has  a  radius  of  0.005 
inches  and  the  center-to-center  spacing  has  been  measured  as  0.050  inches.  These 
values  must  now  be  proportioned  to  give  the  flaw  size  in  a  unit  cube.  The  actual 
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center-to  center  spacing  is  equivalent  to  the  long  diagonal  of  the  RVE  and  the  mea¬ 
sured  flaw  radius  is  proportioned  accordingly.  The  proportionality  Eq.  is 


(radius  of  RVE  flaw) 


n/3~ 


(radius  of  measured  flaw)  (measured  center-to -center  spacing) 


(22) 


For  the  hypothetical  case. 


*J3 


or 


r,  radius  of  RVE  flaw  =  “Q5Q  (.005) 
r  -  0. 173 


(23) 


With  r  now  determined,  and  V  are  easily  calculated.  The  upper  and  lower 
bounds  can  be  determined  if  the  limit  strength  ct  ,  of  the  plastic  matrix  is  known. 

There  are,  of  course,  idealizations  in  the  above  hypothetical  example.  In  real 
materials,  die  flaws  will  not  be  of  a  uniform  size.  Instead,  there  will  be  a  distribu¬ 
tion  of  flaw  sizes.  Which  statistical  measure  of  this  distribution  should  be  used  in 
the  RVE  analysis  can  only  be  determined  from  experiments.  Several  choices  are 
the  mean  size,  the  average  size  or  the  largest  size  in  the  distribution.  In  any  case, 
the  plane  section  to  be  studied  in  a  given  material  will  appear  as  shown  in  Fig.  12(b). 
Further,  the  measurement  of  the  representative  center -to -center  spacing  used  in  the 
hypothetical  example  is  valid  regardless  of  whether  or  not  the  flaw  size  are  uniform 


(23). 

3.3  Limit  Analysis  Applied  to  RVE  of  a  Structure  with  a  Periodic  Array  of  Slits 
The  analysis  in  this  section  is  undertaken  with  the  intent  of  modeling  and  pre¬ 
dicting  failure  in  cast  iron  materials.  Before  discussing  the  RVE  model  chosen 
for  this  case,  it  is  best  to  discuss  the  solidification  process  which  occurs  in  making 
cast  iron.  Figure  13  shows  schematically  and  in  somewhat  simplified  fashion  the 
steps  in  solidification  of  gray  cast  iron  (24).  The  first  step  is  the  formation  of 
graphite  nuclei  throughout  the  molten  metal.  These  nuclei  then  proceed  to  grow 
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arms  in  several  directions.  Next,  the  matrix  metal  between  these  growing  arms  will 
solidify.  Eventually,  the  graphite  arms  from  one  nucleus  approach  arms  from  other 
nuclei  and  the  matrix  material  between  arms  impinge  on  the  matrix  between  arms  of 
adjacen  nuclei.  The  matrix  and  arms  associated  with  one  graphite  nucleus  constitute 
one  eutectic  cell.  The  above  procedure  takes  place  in  three  dimensions;  hence,  the 
eutectic  cell  is  three-dimensional.  In  gray  cast  iron,  the  flakes  are  quite  random 
in  orientation  and  extend  almost  to  the  boundary  of  the  eutectic  cell.  In  nodular  cast 
iron,  the  graphite  flakes  are  condensed  and  form  a  spheroid  in  the  eutectic  cell. 

For  the  purposes  of  analysis,  the  RVE  chosen  is  an  idealized  eutectic  cell  of 
cast  iron.  The  RVE  is  a  cube  of  material  of  unit  dimensions  on  all  sides.  Because 
the  graphite  flakes  in  cast  iron  are  quite  weak,  the  "flakes"  in  the  RVE  will  be 
modeled  as  slits.  These  slits  are  straight  and  they  all  pass  through  the  center  of  the 
RVE.  The  slits  are  oriented  at  various  but  undefined  angles  from  the  center  of  any 
face  of  the  RVE  as  shown  in  Fig.  14.  Figure  14  is  a  sketch  of  1/8  of  the  RVE 
chosen  to  represent  gray  cast  iron.  All  other  segments  of  the  RVE  are  symmetric 
as  indicated  by  the  dashed  extension  lines  in  the  figure.  The  matrix  material  is 
assumed  to  be  the  von  Mises  type  with  tensile  limit  strength  equal  to  a  .  The  slits 
are  assumed  to  have  some  small  but  finite  distance  between  the  slit  faces.  Further, 
the  faces  of  the  slit  will  be  assumed  to  be  very  rough.  The  small  distance  between 
slit  faces  and  the  rough  surfaces  of  the  slit  mean  that  load  or  displacement  will  be 
completely  transmitted  across  the  slit  if  even  a  minute  normal  compressive  load  or 
displacement  occurs  across  the  slit.  Likewise,  if  a  zero  or  normal  tensile  component 
of  load  or  displacement  occurs  across  the  slit,  then  no  load  or  displacement  can  be 
transmitted  across  the  slit.  It  should  be  noted  that  the  sum  total  of  all  flakes  in  the 
RVE  model  of  gray  iron  result  in  a  square  projection  on  any  RVE  face.  For  the 
nodular  iron  case,  the  spherical  nodule  has  a  circular  projected  area  on  any  RVE  face. 
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A  simple  lower  bound  for  the  RVE  can  be  obtained  in  the  same  way  that  the 
lower  bound  was  obtained  for  the  hole  model;  that  is,  assume  that  straight  structural 
members  extend  through  the  matrix  from  one  face  to  its  opposite  face  as  indicated  in 
Fig.  15(a).  Consider  a  plane  cut  through  the  center  of  the  RVE  with  two  opposite  faces 
subject  to  tensile  stress  as  shown  in  Fig.  15(b).  Because  of  the  slits  in  the  structure, 
the  central  portion  of  the  RVE  can  not  carry  tensile  load  and  all  tensile  load  is  car¬ 
ried  by  the  matrix  ligament  at  the  s;de.  Consideration  of  the  total  ligament  area  avail¬ 
able  to  carry  load  perpendicular  to  any  principal  stress  shows  that  the  area  available 
is  just  Am  where  A  is  defined  as  the  total  RVE  face  area  (unity)  minus  the  pro¬ 
jected  area  of  all  slits  onto  the  RVE  face.  If  the  projected  area  of  all  slits  onto  an 
RVE  face  is  a  perfect  square,  then  in  tnaxial  tension  there  can  be  no  uniaxial  stressed 
regions  in  the  structural  members  passed  through  the  matrix.  Hence,  in  triaxial  ten¬ 
sion  there  are  three  biaxial  limit  equations  and  one  triaxial  limit  equation.  These  are 
the  following; 


The  bar  stresses  indicate  that  the  conversion  from  stresses  in  the  structural  members 
to  stresses  on  the  surface  of  the  RVE,  as  discussed  in  the  lower  bound  analysis  of  the 
hole  model,  has  already  been  carried  out  here.  In  three-dimensional  principal  stress 
space,  these  equations  give  a  triaxial  tension  limit  surface  with  a  shape  as  shown  in 
Fig.  16. 
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In  the  case  of  a  circular  projection  of  slits  onto  an  RVE  face  as  in  a  nodular 
cast  iron  model,  there  will  be  uniaxial  stressed  areas  in  triaxial  tension  with  the 
result  that  three  additional  lower  bound  equations  arise  and  govern  the  triaxial  limit 
surface.  These  equations  are  the  following; 

a,  -  A  a 

1  m  my 

Or,  =  A  o 

2  m  my 

and 

o0  -  A  a 
o  m  my 

The  resulting  lower  bound  surface  for  triaxial  tension  in  principal  stress  space  is 
composed  of  three  intersecting  plane  surfaces  as  shown  in  Fig.  17. 

It  is  interesting  to  note  that  if  the  square  projected  shape  of  the  gray  cast  iron 
model  were  slightly  rounded  at  the  corners,  then  uniaxial  stressed  regions  would 
occur  in  triaxial  tension  and  the  resulting  limit  surface  in  principal  stress  space 
would  be  the  same  shape  as  the  nodular  iron  surface,  Fig.  17. 

For  the  stress  situation  involving  two  positive  and  one  negative  principal 
stresses,  the  area  of  the  RVE  surface  acted  on  by  the  positive  stresses  is  still  Am 
but  the  slits  in  the  negative  stress  direction  close  and  transmit  load.  Hence,  the 
area  available  to  carry  load  in  the  negative  direction  is  the  total  RVE  face  area, 
ATot  or  just  unity. 

The  lower  bound  equations  with  6^,  a 2  assumed  positive  and  assumed 
negative  are 


(26a) 

(26b) 
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In  principal  stress  space  with  two  tensile  and  one  compressive  stress,  Eq.  26  appears 
as  shown  in  Fig.  16  or  17. 

In  the  case  of  two  negative  and  one  positive  principal  stresses,  there  are  no 
uniaxially  stressed  regions  in  either  the  cube  or  sphere  RVE.  'The  slit  oriented  per¬ 
pendicular  to  the  tensile  stress  always  opens;  hence,  no  tensile  load  is  transmitted 
through  the  slit  area  of  the  RVE  and  the  area  of  the  RVE  surface  carrying  tensile 


load  is  A  .  The  slits  in  planes  parallel  to  the  compressive  stresses  all  transmit 


load;  hence,  the  total  RVE  face  area  carries  the  compressive  loads.  For  this  case 
then,  the  lower  bound  equations  become  (assuming  and  to  be  the  negative 


stresses), 
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However,  consideration  of  the  various  stress  regions  in  the  assumed  stress 
field  shows  that  the  combinations  of  and  a ^  only  and  and  only  do  not  occur 
anywhere  in  the  assumed  structure  regardless  of  the  shape  of  the  agglomeration  of  slits. 
Hence  in  both  the  gray  iron  and  nodular  iron  models,  the  region  of  the  limit  surface  in 
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which  two  stresses  are  negative  and  one  is  positive  is  governed  by  Eq.  30  for  the 
lower  bound.  This  equation  is  simply  a  von  Mises  ellipse  with  the  center  displaced 
from  the  original  amount  which  depends  on  the  magnitude  of  the  stresses.  This  por¬ 
tion  of  the  limit  surface  is  shown  in  Figs.  16  and  17. 

With  respect  to  Eqs.  24  through  30,  it  should  be  noted  that  cyclic  permutations 
of  the  indices  on  the  stress  terms  give  two  additional  sets  of  equations  in  each  case. 

All  equations  thus  produce  the  complete  three-dimensional  surfaces  as  shown  in 
Figs.  16  and  17. 

The  remaining  stress  condition  to  be  considered  in  the  lower  bound  analysis  is 
the  case  of  triaxial  compression.  In  tnaxial  compression,  all  slits  in  the  structure 
close  up  and  the  area  available  to  carry  load  in  any  of  the  three  directions  is  the  total 
area  of  the  cube  face  or  unity.  Further,  consideration  of  the  stress  state  at  any  point 
in  the  RVE  shows  that  all  points  are  subject  to  triaxial  compressive  stress.  Hence, 
the  only  lower  bound  equation  in  triaxial  compression  is 

7  K  -  32>2  +  <52  -  V2  +  <53 ' =  °m>2  <31) 

This  equation  is  just  the  von  Mises  cylinder  as  shown  in  Figs.  16  and  17. 

Equation  31  completes  the  lower  bound  analysis  for  the  RVE  chosen  to  repre¬ 
sent  gray  and  nodular  cast  irons.  Consideration  of  the  lower  bound  for  the  gray  iron 
case,  shown  in  Fig,  16,  shows  that  there  are  concave  regions  on  the  surface.  While 
it  is  true  that  Drucker  (25)  has  shown  that,  limit  surfaces  must  not  have  concave 
regions,  it  is  noted  that  Fig.  16  is  not  a  true  limit  surface,  but  only  a  lower  bound 
and  hence  may  have  concave  regions. 

Upper  bounds  for  the  proposed  RVE  models  require  that  a  velocity  field  be 
assumed  as  pointed  out  in  (16),  these  velocity  fields  must  be  such  that  RVE  faces 
remain  compatible  with  one  another.  This  restriction  on  any  assumed  velocity  field 
is  necessary  because  RVE's  must  be  able  to  be  stacked  up  to  make  a  larger  structure 
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without  any  gaps  between  them.  If  RVE  faces  were  permitted  to  warp  undei  deforma¬ 
tion,  gaps  would  occur  between  the  warped  face  of  one  RVE  and  that  of  the  next.  This 
is  not  permissible  in  the  total  structure;  hence,  the  restriction  on  the  deformation  of 
the  RVE.  For  the  gray  iron  model,  one  velocity  field  which  satisfies  the  compatibility 
requirements  is  shown  in  Fig.  18(a).  This  ve loci'  '  field  is  shown  on  a  typical  RVE 
cross  section  cut  parallel  to  the  direction  of  an  applied  uniaxial  displacement.  The  uni¬ 
axial  displacement  has  been  assumed  in  the  1  direction  but  the  same  velocity  field  is 
assumed  for  uniaxial  displacement  in  2  or  3  direction.  Further,  portions  of  adjacent 
RVE's  are  shown  in  order  to  indicate  that  compatibility  requirements  are  met. 

The  velocity  field  assumed  in  one  matrix  ligament  of  one  RVE  is  isolated  in 
Fig.  18(b).  Since  the  only  dissipation  in  the  assumed  field  occurs  in  the  von  Mises 
matrix,  the  dissipation  rate  is  simply  two  times  the  yield  in  shear  times  the  shear  dis¬ 
placement  rate  times  the  sheared  area  or 

In  this  equation,  t  is  the  thickness  in  the  plane  of  the  paper  and  .2  is  the  maximum 
length  of  a  diagonal  slit.  For  the  assumed  field,  the  external  work  is  simply  2P^  li. 

The  subscript  1  on  P  indicates  that  the  load  and  displacement  occur  in  the  1  direc¬ 
tion  but  results  are  the  same  for  the  2  and  3  direction.  Setting  the  external  work 
equal  to  the  internal  dissipation  gives  the  upper  bound  limit  load  in  the  1  direction, 
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Substituting  Eq.  33  into  Eq.  32  for  t  gives 


but 
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It  is  noted  that  a  discrepancy  arises  in  the  assumed  plane  strain  velocity  field 
where  the  field  of  one  side  ligament  intersects  that  of  another  side  ligament  at  the  RVE 
edge.  This  discrepancy  will  alter  the  total  internal  dissipation  somewhat  but  as  long  as 
i/42  is  close  to  unity,  the  discrepancy  should  be  small  and  can  be  neglected. 

For  the  extension  of  the  proposed  upper  bound,  Eq.  35,  to  general  three-dimen¬ 
sional  stress  states  involving  at  least  one  tensile  stress,  it  will  be  assumed  that  the 
preceding  velocity  field  occurs  due  to  the  tensile  stress  regardless  of  the  magnitude  of 
any  compressive  principal  stress.  Hence,  Eq.  35  is  a  plane.  Since  the  same  analysis 
holds  for  the  2  and  3  directions  as  for  the  1  direction,  Eq.  35  and  permutations 
thereof  leads  to  three  intersecting  planes  as  shown  in  Fig.  19. 

In  order  to  obtain  a  smaller  upper  bound  in  triaxial  compression  and  in  com¬ 
pression-compression-tension,  it  will  be  assumed  that  the  slits  in  the  structure  close 
and  a  linear  velocity  field  giving  uniform  strains  e^,  e2  and  throughout  the  RVE 
will  occur.  The  dissipation  then  becomes 


ai  el  +  a2  e2  +  a3  c3 


(36) 


where  a^,  a ^  and  a ^  are  stresses  on  the  von  Mises  limit  surface  whert  the  normal 
vector  has  components  proportional  to  e^,  ^  and  The  external  work  is 
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Since  the  strain  rates  are  arbitrary,  the  coefficients  of  each  strain  rate  are  set  equal 
giving 

a.  -  cr.  (38) 

i  l 


However,  a.  must  satisfy  the  von  Mises  limit  ciiterion;  hence, 

1  [<31  •  +  <°2  "  »/  +  <33  '  S1)21  =  amy 

Thus,  the  intersecting  planes  in  Fig.  19  are  cut  off  by  the  von  Mises  cylinder,  Eq.  39, 
as  shown  in  Fig.  19. 

The  upper  bound  shown  in  Fig.  19  has  been  derived  for  the  case  of  the  gray  iron 
model  where  the  slits  define  a  cube  within  the  RVE.  For  the  nodular  iron  model,  the 
slits  define  a  sphere  within  the  RVE.  In  this  case  the  assumptions  made  in  the  gray 
iron  upper  bound  analysis  become  questionable,  particularly  the  assumption  that  the 
change  in  dissipation  due  to  intersecting  deformation  fields  at  the  RVE  edges  is  negli¬ 
gible.  Therefore,  strictly  speaking,  the  lower  and  upper  bounds  for  gray  iron  have 
been  determined  but  only  the  lower  bound  for  nodular  iron  has  been  found.  Since  the 
major  thrust  of  the  experimental  work  in  Chapter  2  and  the  subsequent  discussions  in 
Chapter  4  involve  gray  iron,  the  analysis  of  the  nodular  iron  in  a  theoretical  context 
will  not  be  carried  further. 

Comparison  of  the  upper  and  lower  bound  surfaces  for  the  gray  iron  model, 

Figs.  16  and  19,  reveal  that  in  stress  situations  involving  one  oi  two  negative 
stresses,  the  bounds  are  quite  far  apart.  This  situation  is  observed  more  easily  if 
plane  stress  surfaces  are  compared  as  in  Fig.  20.  Because  of  experimental  evidence 
to  be  pres*  nted  in  the  next  section,  it  is  believed  that  the  discrepancy  in  bounds  is  due 
to  a  poor  velocity  field  in  the  upper  bound  and  that  the  lower  bound  is  generally  closer 
to  tiie  true  limit  load.  Considerable  effort  has  been  put  into  attempting  to  improve  the 


..wims 


upper  bound  in  negative  stress  situations.  This  work  has  met  with  little  success  to 
date  and  remains  an  area  for  future  work. 

Comparison  of  Figs.  16  and  19  also  reveals  regions  of  coincidence.  Three 
points  in  the  tnaxial  tension  octant  have  identical  lower  and  upper  bounds.  This  occurs 
at 
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(40) 


and  permutations  of  this  set  of  equations.  Further,  in  triaxial  compression,  the  lower 
and  upper  bounds  are  identical  von  Mises  cylinders.  Hence,  for  these  stress  states, 
the  exact  limit  condition  has  been  obtained. 

General  three-dimensional  limit  surfaces  are  somewhat  cumbersome  to  work  with 
in  the  comparison  and  analysis  of  other  theories  and  data.  Further,  most  microflaw 
material  failure  theories  available  in  the  literature  are  of  a  two-dimensional,  plane 
stress  nature.  Therefore,  in  the  remainder  of  this  chapter,  the  plane  stress  limit 
surfaces  of  Fig.  20  will  be  used  but  the  existence  of  a  general  stress  state  limit  surface 
is  an  important  aspect  of  the  proposed  theory. 

In  applying  the  proposed  theory  to  actual  materials,  there  are  two  possible  ap¬ 
proaches.  The  first  technique  is  very  much  the  same  as  previously  discussed  for  the 
hole  model.  The  hole  or  flaw  radius  of  the  previous  discussion  is  replaced  by  half  the 
length  of  the  longest  flake  observed  in  a  cross  section  of  the  material.  Interparticle 
center- to -center  spacing  of  the  hole  model  must  be  replaced  by  some  measure  of  the 
eutectic  cell  size  in  the  flake  model.  This  eutectic  cell  size  can  be  obtained  from 
proper  preparation  of  the  polished  surface  of  the  material  or  by  assumption.  With  the 
flaw  size  and  spacing  parameters  determined,  the  real  material  RVE  is  proportioned 
to  the  unit  cube  RVE  by  means  of  Eq.  22  as  before.  Once  this  is  accomplished  and  a 
suitable  estimate  of  the  matrix  material  yield  strength  is  made,  the  upper  and  lower 
bounds  of  the  proposed  theory  are  easily  obtained. 
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Another  and  more  practical  method  of  obtaining  the  parameters  for  the  RVE 
from  an  actual  material  is  to  simply  conduct  uniaxial  tension  and  uniaxial  compression 
tests  on  the  material.  This  information  gives  two  points  on  the  composite  yield  surface. 
The  compression  test  gives  the  strength  of  the  matrix,  crm  via  Eqs.  26,  29  or  31. 

This  point  is  the  same  for  both  upper  and  lower  bounds.  The  tensile  test  gives  an  esti¬ 
mate  of  the  load  carrying  area,  through  Eq.  25.  The  actual  gross  section  failure 

stress  intension  is  between  A  a  and  1.15  A  a  .  Several  tensile  tests  should 

m  my  m  my 

be  conducted  and  the  lowest  value  taken  as  A  a  .  The  highest  value  should  then  be 

m  my  6 

less  than  1.15  A  a  .  Once  a  and  A  are  known,  the  upper  and  lower  bound 
m  my  my  m  ^ 

equations  are  determined. 

This  technique  of  obtaining  the  composite  yield  surface  is  more  practical  because 
metallurgical  structures  can  be  complicated  in  microflawed  materials,  and  determining 
the  microscale  parameters  via  photomicrographs  can  be  difficult. 

3.4  Comparison  of  Proposed  Failure  Theory  with  Currently  Available  Theories 

In  the  past,  numerous  theories  for  failure  of  microflawed  materials  have  been 
developed.  Perhaps  the  simplest  of  these  failure  theories  was  that  of  Fisher  (26). 

This  investigator  hypothesized  that  plastic  action  occurred  at  the  tip  of  graphite  flakes 
in  gray  cast  iron  at  very  small  loads  and  that  fracture  followed  soon  after.  The  dis  ¬ 
tortion  energy  or  von  Mises  criterion  for  yielding  was  thought  to  best  characterize 
this  plastic  action  and  it  was  hypothesized  that  fracture  could  also  be  characterized 
by  this  criterion.  Fisher  further  assumed  that  in  tension,  microflaws  opened  and 


caused  stress  concentrations  in  the  matrix  while  in  compression,  the  microflaws 
tended  to  close  and  transmit  load.  Fisher's  aim  was  to  predict  the  failure  of  gray 
cast  iron;  hence,  his  model  was  slanted  toward  flake  type  microflaws  rather  than 
spherical  types.  Nevertheless,  the  equations  he  developed  are  general  and  can  be 
applied  to  either  case. 
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Coffin  (27)  considered  Fisher's  theory  to  provide  a  lower  bound  to  the  fracture 
stress  of  gray  cast  iron  and  hypothesized  that  actual  failure  stresses  should  be  higher. 
Therefore,  Coffin  developed  a  failure  criterion  based  on  the  stress  concentration  at  the 
tip  of  graphite  flakes  in  gray  cast  iron.  According  to  Coffin's  theory,  there  is  a  resid¬ 
ual  stress  at  the  tip  of  graphite  flakes  which  arises  due  to  the  different  thermal  expan¬ 
sion  characteristics  of  the  graphite  and  matrix. 

In  general.  Coffin's  equations  for  failure  are  more  complex  than  Fisher's  and 
require  more  constants.  Further,  there  seems  to  be  little  support  for  the  magnitude 
of  the  residual  stress  term  that  is  required  to  fit  experimental  data.  Coffin  found  that 
a  residual  stress  of  30  ksi  was  required  to  fit  his  equation  to  experimental  data. 

Clough  and  Shank  (28)  later  showed  that  this  magnitude  for  the  residual  stress  is  very 
unlikely. 

In  a  study  of  masonry  failure.  Coulomb  (29)  found  that  failure  occurred  when  the 
shear  stress  on  some  particular  plane  reached  a  critical  value.  The  particular  plane 
of  importance  in  his  equation  is  a  material  dependent  property.  Paul  (30)  later  showed 
that  the  Coulomb  theory  incorrectly  predicted  the  orientation  of  the  fracture  surface  in 
biaxial  tension  and  in  combined  tension-  compression.  This  investigator  proposed  that 
fracture  will  occur  when  either  the  Coulomb  criteria  or  the  maximum  tensile  strength 
of  the  material  is  reached  and  whichever  is  reached  first  is  when  failure  occurs.  Thus, 
the  Coulomb  theory  is  combined  with  a  maximum  principal  stress  criterion  to  give  a 
fracture  theory  which  Paul  found  can  accurately  predict  the  fracture  plane  and  the  stress 
at  fracture. 

Mann  (31)  revised  the  Coulomb  theory  by  proposing  a  quadratic  relation  mvolv 
mg  the  critical  shear  and  normal  stresses  on  some  particular  plane  rather  than  a 
linear  relation  as  originally  proposed  by  Coulomb.  Marin's  theory  generally  fits  data 
better  than  Coulomb's  but  suffers  from  the  fact  that  three  separate  material  dependent 
properties  are  required. 


Babel  and  Sines  (32)  developed  a  fracture  criterion  for  a  microflawed  material 
which  involves  the  aspect  ratio  of  an  assumed  two-dimensional  elliptical  flaw.  The 
basis  of  this  criterion  developed  by  Babel  and  Sines  was  stress  concentration  equation 
at  various  locations  around  an  elliptical  flaw  which  were  developed  by  Inglis  (33).  No 
attempt  was  made  in  this  theory  to  incorporate  a  parameter  based  on  the  spacing  of 
the  flaws.  The  theory  of  Babel  and  Sines  predicts  that  for  circular  flaws  in  a  brittle 
matrix  the  compressive  failure  stress  will  be  three  times  the  tensile  stress  while  for 
long,  thin  elliptical  shape  flaws  the  compressive  to  tensile  failure  stress  ratio  ap¬ 
proaches  seven. 

Figure  21  shows  the  above  mentioned  theories  assuming  a  tensile  stress  of 
33  ksi  and  a  compressive  stress  of  100  ksi.  All  theories  except  that  of  Babel  and 
Sines  are  identical  in  biaxial  tension  and  agree  well  with  experimental  data.  In 
tension-compression.  Coffin  and  Fisher's  theories  are  coincident  while  Paul’s 
theory  is  higher  in  the  tension-compression  region.  In  biaxial  compression, 

Fisher  s  theory  in  plane  stress  is  just  the  von  Mises  ellipse;  Coulomb  and  Marin’s 
theories  are  straight  lines.  The  theory  of  Babel  and  Sines  does  not  include  biaxial 
compression  estimates,  so  this  theory  1S  not  extended  into  this  quadrant. 

Several  comments  on  the  currently  available  failure  theories  for  microflawed 
materials  are  in  order.  Most  of  the  theories  discussed  involve  numerous  "material 
dependent  properties  1  or  constants.  The  theories  of  Babel  and  Sines  and  Fisher  have 
the  least  number  of  constants  with  one.  However,  the  Babel  and  Sines  theory  does 
not  include  biaxial  compression  predictions.  All  of  the  currently  available  theories 
have  been  developed  from  two  dimensional  models  and  cannot,  in  general,  be  ex¬ 
tended  m  a  rigorous  manner  to  three-dimensional  stress  conditions.  Further,  none 
of  the  theories  discussed  involve  a  flaw  spacing  parameter  explicitly  and,  hence- 
would  predict  the  same  material  behavior  for  all  microflawed  materials  regardless 
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or  flaw  spacing.  Finally,  only  the  theory  of  Babel  and  Sines  includes  a  factor  in¬ 
volving  the  shape  of  the  microflaw. 

Comparison  of  the  proposed  theory  for  the  slit  model  with  the  theories  just 
discussed  shows  that  the  lower  bound  equation  is  identical  to  Fisher's  theory.  This 
supports  Coffin's  view  that  Fisher's  theory  is  a  lower  bound  to  failure.  Further 
comparison  of  the  various  theories  can  be  made  by  reconsidering  Fig.  21  with  the 
proposed  theory  upper  and  lower  bounds  drawn  in.  In  biaxial  tension  and  combined 
tension -compression,  all  existing  theories  are  encompassed  by  the  proposed  bounding 
equations.  In  biaxial  compression,  the  proposed  upper  and  lower  bounds  are  coinci¬ 
dent  with  each  other  and  with  the  von  Mises  ellipse  of  the  matrix  material  alone. 

The  essential  differences  between  the  proposed  theory  and  those  currently 
available  are  that  the  proposed  theory  incorporates  measures  of  flaw  size,  shape 
and  spacing  whereas  none  of  the  other  theories  include  all  these  factors  explicitly. 
Further,  the  proposed  theory  involves  a  general  three-dimensional  model  and  has 
been  extended  in  a  rigorous  manner  to  include  general  triaxial  stress  conditions. 

_ Comparison  of  Proposed  Theory  with  Experimental  Data 

Comparison  of  the  proposed  theory  for  the  hole  model  with  experimental  data 
on  this  type  of  material  is  not  practical  because  there  are  no  materials  available 
which  approximate  the  hole  model  throughout  a  given  structure  in  a  homogeneous 
and  statistically  isotropic  manner.  Sintered  powder  metallurgy  products  come  the 
closest  to  the  model  but  can  not  be  used  for  two  reasons.  First,  the  objective  of 
the  manufacturer  of  sintered  powder  products  is  to  decrease  the  volume  fraction  of 
voids.  Hence,  the  volume  fraction  of  voids  in  commercially  available  material  is 
quite  small,  often  less  than  2%.  Secondly,  because  of  the  sintered  powder  process, 
suifaces  of  manufactured  parts  are  very  dense  and  porosity  is  concentrated  in  the 
center.  The  hole  model  of  the  proposed  theory  could  be  used  to  model  the  various 


layers  is  a  sintered  powder  product  by  changing  the  size  and  spacing  of  the  RVE  for 
various  layers.  However,  the  analysis  of  the  total  layered  structure  then  becomes  an 
important  but  solvable  problem  in  its  own  right.  This  type  of  analysis  is  considered 
to  be  beyond  the  scope  of  the  current  study  and  further  attention  to  the  hole  model  is 
not  warranted  here. 

For  the  case  of  the  slit  model,  the  model  was  set  up  specifically  to  model  cast 
iron  and,  therefore,  comparison  of  the  proposed  theory  with  experimental  data  on 
cast  irons  should  be  conducted.  Two  types  of  cast  irons  will  be  compared  with  the 
proposed  theory.  These  are  gray  and  nodular  cast  irons.  Figure  22  shows  experi¬ 
mental  data  on  fracture  of  gray  iron  tubes  under  various  ratios  of  applied  stress  after 
Cornet  and  Grassi  (35).  Also  shown  in  this  figure  is  che  Fisher  theory,  the  modified 
Coulomb-Mohr  theory  of  Paul,  and  the  upper  and  lower  bound  of  the  proposed  theory. 
The  proposed  theory  constants  were  determined  from  the  tensile  and  compressive 
fracture  stress  using  25  ksi  as  the  tensile  failure  stress  and  90  ksi  as  the  ultimate 
stress  of  the  matrix  alone.  The  Fisher  theory  gives  a  good  estimate  of  the  failure 
stress  in  biaxial  tension  but  underestimates  the  failure  in  tension-compression.  The 
modified  Coulomb  theory  of  Paul  fits  the  data  slightly  better  than  Fisher's  theory. 

The  proposed  theory  gives  upper  and  lower  bounds  which  are  15%  apart  in  biaxial 
tension  but  as  much  as  501/,  apart  in  tension -compression. 

Figure  23  shows  fracture  data  on  gray  cast  iron  from  Coffin  (27)  along  with 
Coffin's  theoretical  fracture  criterion.  Also  shown  in  this  figure  are  the  Fisher 
theory  and  the  upper  and  lower  bounds  of  the  proposed  theory.  Coffin  and  Fisher's 
theory  are  coincident  in  biaxial  tension  and  in  tension  and  compression.  I  isher  s 
theory  in  plane  stress,  and  the  upper  and  lower  bounds  of  the  proposed  theory  are 
coincident  in  biaxial  compression.  The  lower  bound  of  the  proposed  theory  was 
based  on  a  matrix  ultimate  strength  of  110  ksi  and  an  Am  of  0.254.  These  values 
were  chosen  to  fit  the  experimental  data  in  uniaxial  tension  and  compression.  The 
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lower  bound  fits  the  data  fairly  well  but  not  as  well  as  the  Coffin  or  Fisher  theory. 

The  upper  bound  is  far  from  the  data  in  tension-compression.  In  biaxial  compression, 
the  upper  and  lower  bounds  are  coincident  and  they  do  not  appear  to  fit  the  experimental 
data.  This  discrepancy  is  attributed  to  difficulty  in  conducting  a  biaxial  compression 
test  (34). 

In  a  study  involving  experimental  data  on  gray  cast  iron,  Mair  (34)  found  that 
Fisher's  theory  was  the  easiest  to  use  and  gave  the  best  overall  agreement  with  the 
data.  Figure  24  shows  Mair's  comparison  (34)  in  addition  to  the  bounds  of  the  pro¬ 
posed  theory.  Mair  discussed  the  importance  of  the  third  principal  stress  in  biaxial 
compression  and  Fig.  24  shows  Fisher's  theory  in  both  plane  stress  and  in  triaxial 
stress  states.  The  consideration  of  the  third  principal  stress  in  Fisher’s  theory  fits 
the  experimental  results  much  better  than  the  plane  stress  theory.  This  fact  points 
out  the  importance  of  the  third  stress  in  biaxial  compression.  It  will  again  be  noted 
that  the  upper  bound  of  the  proposed  theory  overestimates  the  actual  failure  in  the 
tension -  compression  quadrant. 

Figure  25  shows  experimental  results  of  Clough  and  Shank  (28)  along  with 
Fisher's  theory  and  the  bounds  of  the  proposed  theory.  In  developing  Fisher's  theory 
and  the  proposed  theory  lor  this  data  a  value  of  140  ksi  was  used  for  uniaxial  com  ¬ 
pression  and  40  ksi  for  umaxial  tension  failure.  Again,  Fisher's  theory  fits  the  data 
fairly  well,  and  the  upper  bound  of  the  proposed  theory  overestimates  the  tension- 
compression  failure. 

Although  the  proposed  theory  gives  a  rigorous  lower  bound  but  no  upper  bound 
for  the  nodular  iron  model,  comparison  of  the  lower  bound  with  experimental  data  can 
provide  interesting  information. 

Figure  26  shows  experimental  results  of  Cornet  and  Grassi  (36)  for  nodular 
iron  along  with  the  bound  of  the  proposed  theory.  The  matrix  ultimate  strength  was 
taken  as  75  ksi  and  the  load  carrying  area  was  taken  as  0.  748  square  inches.  These 
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values  were  chosen  to  fit  the  experimental  data.  The  lower  bound  of  the  proposed 
theory  fit  the  data  well. 

Figure  26  also  shows  experimental  data  from  Clough  and  Shank  (37)  for  nodular 
iron.  Again,  the  constants  for  the  lower  bound  of  the  proposed  solution  were  chosen  to 
fit  the  experimental  data.  The  lower  bound  for  the  proposed  theory  again  fit  the  data 
well. 

All  of  the  examples  used  for  comparison  thus  far  have  involved  fitting  the  pro¬ 
posed  theory  to  experimental  data  by  using  the  macroscale  uniaxial  tension  and  com¬ 
pression  data.  This  is  because  the  proposed  theory  requires  an  accurate  determina¬ 
tion  of  the  matrix  strength  as  well  as  the  spacing  and  size  of  flaws  in  the  matrix.  In 
the  investigations  previously  cited,  one  or  more  of  these  parameters  could  not  be 
accurately  determined.  However,  in  a  recent  study  on  nodular  cast  iron,  Testin  (38) 
measured  the  strength  of  the  matrix  material  as  well  as  giving  photomicrographs  of 
the  structure  of  three  nodular  cast  irons.  The  matrix  strength  of  the  irons  was  171  ksi. 
Three  different  nodule  sizes  were  studied.  The  largest  nodule  size  and  average  center- 
to-center  spacing,  as  measured  from  the  photomicrographs  of  Testin,  are  shown  in 
Table  Id  The  nodule  size  and  spacing  is  proportioned  to  a  unit  cell  size  for  the  model 
as  discussed  in  an  earlier  section.  The  equivalent  diameter  of  the  nodule  in  the  unit 
cell  is  listed  in  Table  10  along  with  the  load  carrying  area.  The  lower  bound  is  the 
matrix  ultimate  stress  multiplied  by  Am>  The  measured  failure  stresses  for  each  of 
the  nodule  sizes  tested  are  listed  in  Table  10,  The  actual  failure  stress  is  quite  close 
to  the  lower  bound  prediction  for  the  small  and  medium  nodule  sizes,  but  the  actual 
stress  is  much  higher  than  the  lower  bound  for  the  large  nodule  size.  This  discrep¬ 
ancy  in  the  last  case  is  partially  attributable  to  the  fact  that  the  photomicrograph  from 
which  the  nodule  size  and  spacing  were  measured  contained  very  few  nodules.  Hence, 
it  was  not  possible  to  determine  a  statistically  significant  nodule  spacing. 
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The  above  correlations  of  the  proposed  theory  with  experimental  data  indicate 
that  the  lower  bound  is  certainly  a  valid  criteria  to  use  for  predicting  lower  bounds 
on  failure  stress  under  various  biaxial  stress  conditions  in  cases  where  external 
notches  are  not  present.  In  the  gray  iron  model,  the  upper  bound  generally  over¬ 
estimates  the  failure  stress,  especially  in  the  biaxial  tension -compression  case. 

Some  currently  available  theories  predict  failure  stresses  better  than  the  proposed 
theory,  but  these  theories  do  not  permit  prediction  of  failure  stresses  from  micro¬ 
scale  quantities  alone  as  does  the  proposed  theory.  This  ability  of  the  proposed 
theory  was  illustrated  in  the  comparison  with  Testin’ s  data  on  nodular  iron. 

Comparison  of  the  proposed  theory  with  investigations  aimed  at  determining 
the  effect  of  casting  variables  on  failure  is  also  of  interest.  For  example,  the  theory 
developed  here  predicts  that  the  strength  of  cast  iron  depends  on  the  size  of  the  graphite 
flakes  in  relation  to  the  eutectic  cell.  If  this  relation  remains  the  same,  the  strength 
of  the  iron  is  constant  regardless  of  the  absolute  size  of  the  graphite.  However,  ex¬ 
perimental  evidence  (39)  shows  that  cast  iron  strength  increases  as  graphite  absolute 
size  decreases.  This  apparent  anomaly  can  be  explained  by  considering  more  care¬ 
fully  the  solidification  process  in  production  of  fine  graphite.  Fine  graphite  is  pro¬ 
duced  by  cooling  molten  iron  very  quickly.  However,  this  process  leads  to  large 
austenite  dendiites.  Solidification  of  the  fine  graphite  occurs  after  those  dendrites 
are  formed  and  between  the  dendrite  arms.  Hence,  the  theory  developed  herein  is 
not  applicable  to  the  total  structure  in  this  case  but  only  to  the  fine  graphite  area. 

The  theory  developed  must  then  be  restricted  to  usage  in  cases  where  graphite  is 
large,  completely  random  in  spacing  and  the  austenite  dendrites  are  small.  This 
implies  use  only  in  cases  of  slowly  cooled  cast  iron. 

The  existence  of  von  Mises  behavior  in  cast  iron  under  triaxial  compression 
loading  is  predicted  by  the  proposed  theory  and  has  been  studied  to  some  extent  by 
Brandes  (40).  This  investigator  subjected  nodular  cast  iron  to  a  hydrostatic  com- 
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pressive  stress  of  -136  kg/mm  then  applied  tensile  stress  to  one  axis  until  failure 
occurred.  Application  of  the  three-dimensional  form  of  the  von  Mises  theory  with 
o  and  a„  equal  to  -136  kg/mm  and  the  materials  tensile  strength  equal  to  Brandes 

JL  ^ 

2 

uniaxial  tensile  strength  (48.4  kg/mm  )  gives 

■*> 

(48. 4)2  =  i  [(-136  -  a3)2  +  (Og  +  136)2] 

or  (46) 

2 

a3  -  -87  kg/mm 

2 

Brandes  found  that  the  Og  stress  was  -84  kg/mm  for  failure.  Although  this  is  only 
one  test,  it  does  indicate  that  the  predicted  behavior  from  the  proposed  model  in  tri- 
axial  compression  is  reasonable. 

Now  that  the  limit  surface  for  the  microflawed  material  without  macroscale 
notches  has  been  obtained  and  its  validity  verified  with  experimental  results,  the 
case  of  macroscale  notches  can  be  studied.  This  extension  of  the  proposed  theory 
comprises  Chapter  4. 
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4.  MACROSCALE  LIMIT  ANALYSIS  OF  A  M1CROFLAWED  MATERIAL 
4. 1  Introduction 

Now  that  the  limit  behavior  of  a  micioscale  representative  element  has  been 
obtained,  the  limit  surface  can  be  used  in  macroscale  analyses.  As  previously 
stated,  since  the  modeled  element  is  representative,  any  structure  can  be  built 
up  by  stacking  elements;  hence,  the  RVE  limit  surface  can  be  applied  directly  to 
the  macroscale  problem.  The  macroscale  limit  analysis  is  conducted  by  using  the 
microscale  RVE  limit  condition  developed  herein  as  the  material  continuum  limit 
condition.  This  extension  is  valid  under  the  conditions  set  by  Majumdar  and 
McLaughlin  (16)  and  discussed  in  Chapter  3. 

Three  types  of  notched  specimens  will  be  studied.  These  are  the  double  edge 
cracked  plate,  the  circumferential  notched  bar,  and  the  notched  three  point  bend  bar. 

The  analyses  will  employ  the  RVE  limit  surfaces  of  the  flake  model  obtained  in 
Chapter  3.  Specifically,  the  limit  surfaces  of  the  gray  iron  model.  Fig.  16  and  19 
will  be  used.  In  order  to  reduce  the  number  of  calculations  in  this  section,  only  an 
estimated  single  Emit  surface  will  be  used  here.  This  estimated  surface  is  shown 
in  Fig.  27  and  represents  a  limit  surface  between  the  established  lower  and  upper 
bounds  of  Chapter  3. 

The  results  of  the  macroscale  analyses  of  the  various  specimens  will  be 
compared  with  the  experimental  data  from  the  current  study  presented  in  Chapter  2. 


The  configuration  of  the  specimen  to  be  studied  here  is  shown  in  Fig.  6.  The 
plate  width  is  W,  the  thickness  is  t,  and  the  crack  depth  is  "a"  from  each  side  of 
the  plate.  A  simple  lower  bound  is  obtained  by  considering  a  structural  member  to 
extend  through  the  center  of  the  plate  as  shown  in  Fig.  28.  The  lower  limit  load  of 
the  structure  is  the  limit  load  of  this  member  or 
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PLB  =  °oc  (W  "  2a>  1  (42) 

where  is  the  lower  bound  limit  load,  and  is  the  limit  stress  of  the  micro- 
flawed  material  in  uniaxial  tension.  An  upper  bound  for  the  limit  load  can  be  obtained 
by  considering  a  normal  displacement  velocity  field  as  shown  in  Fig.  28.  This  assumed 
field  requires  no  strain  in  the  other  two  directions  of  the  plate.  Referring  to  Fig.  27, 
the  uniaxial  strain  assumed  implies  rate  vector  is  normal  to  point  L.  Hence,  the 
internal  energy  dissipated  by  the  assumed  displacement  is 

/  aij  *ij  dV  =  2* 15  CT0C  “  *<W  *  2a>  (43) 

V 

The  external  work  is  the  external  load,  Pyg,  multiplied  by  the  external  displacement  or 

PUB  “  (44) 

Setting  the  external  work  equal  to  the  internal  work  gives 

PUB  =  1,15  °oc  t(W  "  2a)  (45) 

In  this  case,  the  upper  and  lower  bounds  are  within  15%  of  each  other. 

In  this  expression  developed,  oqc  is  the  limit  stress  of  the  RVE  in  uniaxial 
tension.  Using  a  oqc  of  42.5  ksi,  which  is  the  ultimate  stress  of  test  P-M-2  from 
Table  4,  and  the  specimen  dimensions  of  DEC  3  and  4  from  Table  8,  the  two  bounds  are 

Pgg  -  42.  5  (1)  (1. 12)  =  47.6  kips  (46) 

and 

PUB  =  1* 15  (PLB^  =  54'  76  kips  (47) 

From  Table  8,  the  fracture  loads  for  these  two  specimens  were  43.9  and  47  kips. 

Thus,  the  actual  failure  loads  lie  below,  but  close  to,  the  predicted  bounds  of  limit 
analysis. 


Returning  to  the  linear  elastic  fracture  mechanics  techniques  and  using  K^c  = 
18.  8  ksiN/Im,  a  failure  load  prediction  of  25.3  kips  is  obtained.  This  value  under¬ 
estimates  the  failure  load  of  the  DEC  specimens.  Thus,  it  appears  that  gray  cast 
iron  with  macroscale  notches  may  be  too  ductile  to  utilize  LEFM  techniques  but  not 
ductile  enough  to  enable  the  use  of  limit  analysis.  Nevertheless,  LEFM  provides  a 
lower  bound  to  failure  while  limit  analysis  gives  an  upper  bound. 

4.3  Macroscale  Analysis  of  a  Circumferential  Notched  Round  Bar 

The  analysis  for  this  case  is  almost  the  same  as  for  the  double  edge  cracked 
plate;  hence,  the  analysis  will  not  be  repeated.  The  comparison  with  the  data  of  the 
current  study  can  be  accomplished  with  one  minor  change  in  the  previous  equation. 

9 

Instead  of  B(W  -  2a)  in  the  bound  equation,  the  expression  7t/D  where  D'  is  the 
diameter  at  the  notch  must  be  substituted  in  the  double  edge  cracked  plate  equation. 
From  Table  7,  the  three  diameters  tested  were  0.354,  0.  707  and  1.59  inches. 

Figure  29  is  a  plot  of  the  specimen  area  at  the  notch  versus  the  failure  load.  The 
data  points  are  for  the  specimens  tested,  the  solid  lines  are  on  LEFM  methods. 

In  the  low  load,  small  cross  sectional  area  region  of  the  plot,  the  data  and  the  two 
theoretical  techniques  are  quite  close  together  but  as  specimen  size  increases,  the 
LEFM  technique  provides  a  lower  bound  to  failure  and  the  proposed  limit  analysis 
approach  gives  an  upper  bound.  It  is  noted  that  the  actual  failure  loads  follow  the 
general  trend  of  the  fracture  mechanics  prediction  rather  than  the  limit  analysis 
prediction.  Thus,  for  tensile  loading  of  notched  bars  of  gi  jy  cast  iron,  LEFM 
techniques  provide  a  close  but  conservative  estimate  of  failure  loads. 

4.4  Macroscale  Analysis  of  a  Notched  Three  Point  Bend  Bar 

In  conducting  the  lower  bound  analysis  of  three  point  bend  bars,  the  3PB  con¬ 
figuration  can  be  replaced  by  a  pure  bending  configuration  with  moment,  M  equal  to 
PS/4  where  P  is  the  load  from  3PB  bar  and  S  is  the  span  length  in  the  3PB  specimen. 
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This  external  moment  must  be  resisted  by  the  internal  moment  set  up  by  the  fully 
plastic  tensile  and  compressive  stresses  as  shown  in  Fig.  30.  The  tensile  and  com¬ 
pressive  loads  must  be  equal  but  since  the  tensile  and  compressive  limit  stresses  are 
not  equal,  the  areas  of  tensile  load  and  compressive  load  can  not  be  equal.  Hence,  as 
shown  in  Fig.  30,  the  tensile  load  is 

PT  =  aoc  (c)  C  <48) 

and  the  compressive  load,  is 

PC  =  ac  (B  '  C)  t  (49) 


where  C  is  depth  of  beam  subject  to  tensile  stress,  B  is  the  total  depth  of  the  beam 
at  the  notch  and  t  is  the  thickness  of  the  beam.  Setting  PT  =  Pc  and  letting  R  = 

V^Comp  ®lves 


c  = 


B 

(l  +  R) 


(50) 


The  internal  moment  is 


M  = 


(51) 


Setting  the  internal  moment  equal  to  the  external  moment,  gives  the  lower  bound  limit 
load 


PLB  S  jl  +  Rj  1  (52> 

An  upper  bound  is  somewhat  more  difficult  to  obtain  when  the  material  has  dif¬ 
ferent  tensile  and  compressive  limit  loads.  For  the  upper  bound,  a  collapse  mode 
must  be  developed.  The  assumed  collapse  mode  is  shown  in  Fig.  30.  The  energy 
dissipated  in  the  deforming  triangular  region  in  the  center  of  the  bar  can  be  easily 


calculated  from 


p  wjivuvr»>?  l  j  i_h A-_pi  y.41*  V 
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a, .  e. .  dV 
ij  U 


The  strain  in  the  deforming  region  is  uniaxial  tension.  Hence, 


cr. .  e  .  =  cr  t  (1.  15) 
1]  1]  oc  1 


but. 


el  = 


2AB 


(53) 


(54) 


(55) 


for  a  unit  base  of  the  triangular  deforming  region.  A  is  the  maximum  beam  velocity. 
The  volume  integral  (53)  in  this  case  reduces  to 


AB2  t  cr 


oc 


(2.30) 


(56) 


The  exter  nal  work  done  is  P  *  A  so  that  the  upper  bound  load  PUB  is 


B2  t  a 


oc 


UB 


(2.30) 


(57) 


The  ratio  of  the  lower  to  upper  bounds  is 

sjf  ■  {ttr}  <°- 87)  <58) 

The  next  step  is  to  compare  the  limit  load  prediction  with  experimental  data 
from  this  study.  For  the  gray  cast  iron  under  study,  R  =  1/3  and  oqc  is  42.5  ksi. 

B  is  the  depth  of  the  beam  at  the  notch.  Using  S  -  12.  89  from  Table  5,  and  calcu¬ 
lating  the  limit  load  for  various  B  values,  a  trend  is  evident.  Figure  31  is  a  plot  of 
crack  length  a,  which  is  simply  (3.22  -  B)  versus  the  limit  load  stress.  The  solid 
lines  are  the  lower  bound  using  oqc  of  42.  5  ksi  and  the  upper  bound  using  1. 53  times 
the  lower  bound.  The  actual  data  points  from  all  tests  listed  in  Table  5  are  indicated 
as  well  as  the  predictions  of  LEFM  based  on  a  fracture  toughness  of  18.  8  ksiN/In. 
Again,  the  LEFM  predictions  fall  below  the  data  while  the  limit  analysis  approach 
provides  an  upper  bound. 
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4.5  Conclusions  of  Macroscale  Analysis 

From  the  comparison  of  experimental  data  with  the  LEFM  techniques  and 
limit  analysis  predictions,  it  can  be  concluded  the  LEFM  method  should  provide 
lower  bound  predictions  of  failure  in  notched  gray  iron  parts  while  limit  analysis 
should  always  provide  an  upper  bound.  The  implication  of  this  conclusion  is  that 
gray  cast  iron,  in  a  macroscale  notched  situation,  is  too  ductile  to  validate  strict 
use  of  LEFM  but  not  ductile  enough  to  permit  the  strict  application  of  limit  analysis. 

McLaughlin  (22),  has  shown  that  limit  analysis  will  always  provide  an  upper 
bound  to  the  failure  load  of  a  material  which  lacks  sufficient  ductility  for  strict  appli¬ 
cation  of  limit  theory.  Indeed,  this  view  is  supported  by  the  results  of  the  current  study. 

In  order  to  obtain  a  better  estimate  of  the  actual  failure  load  in  cast  iron  parts, 
two  alternatives  appear.  The  first  is  to  modify  LEFM  methods  while  the  second  is 
to  modify  limit  analysis  techniques  to  account  for  limited  ductility.  The  first  method 
is  under  study  by  others  and  is  subject  to  the  objections  raised  in  Chapter  2.  The 
second  alternative  may  prove  useful  and  a  simple  technique  can  be  used  to  enable 
better  prediction  of  the  3PB  test  data. 

Although  it  must  be  conceded  that  limit  analysis  is  not  strictly  applicable  to 
gray  cast  iron  in  bending,  the  remaining  beam  depth,  B  can  be  artificially  reduced 
by  some  amount,  5  to  account  for  the  lack  of  tensile  ductility  and  subsequent  pre- 
catastrophic  crack  extension.  This  6  value  will  be  material  dependent  but  should  be 
geometry  independent.  In  order  to  determine  a  6  for  the  case  at  hand,  a  crack  length 
of  0.5  inches  (B  =  2.  72)  will  be  used.  Further,  assume  that  the  desired  lower  bound 
limit  load  is  16.  5  kips  for  the  case  of  a  -  0,  5.  Using  Eq.  52  and  solving  for  (B  -  6) 
where  (B  -  6)  takes  the  place  of  B,  gives 


(B  *  6)  =  l.  824 


(59) 
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but  B  is  originally  2.72;  therefore, 

6-0. 8957  (60) 

This  is  the  amount  of  crack  extension  which  should  be  added  to  the  original  crack 
length,  a,  in  order  to  have  the  limit  analysis  solution  match  the  actual  data.  If  this 
6  is  a  material  property,  one  should  be  able  to  take  all  of  the  3PB  tests,  add  6  to 
the  actual  crack  lengths  and  predict  failure  using  (B  -  6)  in  Eq.  52.  This  exercise 
has  been  carried  out  and  Fig.  31  shows  the  results.  Note  that  the  actual  data  fall 
either  within  the  limits  or  reasonably  so. 

The  above  modification  of  limit  analysis  remains  to  be  proven  more  exten¬ 
sively.  There  are  a  number  of  important  questions  which  must  be  answered.  For 
example,  is  the  6  factor  geometry  independent  or  will  it  change  with  beam  depth  or 
thickness?  An  important  theoretical  question  is,  ‘Is  the  modified  remaining  ligament 
fully  plastic  at  failure?  ‘  If  not,  use  of  limit  analysis  is  on  poor  footing.  The  con¬ 
clusion  and  use  of  the  above  operational  procedure  must  be  viewed  as  preliminary  in 
nature  at  this  time.  Nevertheless,  the  procedure  is  interesting  and  certainly  can  be 
used  to  predict  the  results  in  the  current  case. 
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5.  SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE  STUDY 


5. 1  Summary  and  Conclusions 

It  has  been  shown,  using  experimental  data  gathered  from  various  specimen 
geometries,  that  Unear  elastic  fracture  mechanics  cannot  be  used  to  predict  failure 
in  gray  cast  iron  independent  of  geometry,  h  may  be  possible  to  predict  precata- 
strophic  crack  extension  or  provide  a  lower  bound  estimate  of  failure  by  using  plane 
strain  fracture  toughness,  but  prediction  of  accurate  failure  loads  is  not  possible. 

It  has  also  been  shown  that  specimen  and  test  criteria  currently  used  in  plane  strain 
fracture  toughness  testing  lead  to  discrepancies  in  the  case  of  cast  iron. 

in  an  effort  to  provide  a  failure  prediction  technique  for  cast  iron,  a  study  of 
limit  analysis  techniques  applied  to  general  microflawed  materials  with  ductile 
matrices  was  undertaken.  The  method  employed  was  to  develop  general  limit  sur¬ 
faces  for  microscale  representative  elements  of  the  material.  These  limit  criteria 
were  then  used  in  macroscale  analysis.  The  proposed  theory  was  compared  with 
currently  available  theories  for  microflawed  material  as  well  as  experimental  data 
on  smooth  specimen  tests  of  gray  and  nodular  cast  irons.  Most  theories  and  exper¬ 
imental  data  for  smooth  specimen  cases  fell  within  the  bounds  of  the  proposed  theory. 
Generally,  those  experimental  results  falling  outside  the  proposed  theory  are  attrib¬ 
utable  to  experimental  difficulty  in  combined  stress  testing  of  smooth  specimens. 

The  proposed  theory  differs  from  others  in  two  major  respects.  First,  the  proposed 
theory  includes  measures  of  microflaw  size,  shape,  spacing,  and  matrix  strength 
properties  explicitly,  whereas  other  theories  include  these  factors  implicitly  through 
macroscale  measures  such  as  tensile  and  compressive  strength.  Also,  the  proposed 
theory  has  been  developed  for  general  three  dimensional  stress  states  where  most 
other  theories  have  been  applied  to  plane  stress  only. 
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The  limit  surfaces  obtained  from  the  microscale  RVE  model  of  gray  cast  iron 
has  been  extended  to  the  macroscale,  and  three  macroscale  fracture  toughness  speci¬ 
mens  have  been  analyzed  using  limit  analysis.  The  results  of  this  analysis  and  com¬ 
parison  with  data  indicate  that  limit  analysis  techniques  provide  an  upper  bound  to  the 
failure  load  of  gray  cast  iron  parts  with  notches.  In  the  case  of  three  point  bend  bars, 
the  limit  analysis  so’uuon  may  be  modified  by  a  material  dependent  constant  to  permit 
prediction  of  failure  loads  for  this  geometry. 

The  overall  conclusions  of  this  study  are  as  follows; 

(1)  Linear  elastic  fracture  mechanics  methods  will  give  lower  bound 
predictions  for  gray  cast  iron  failure  loads  in  cases  where  macro¬ 
scale  notches  are  present,  However,  gray  cast  iron  is  too  ductile 
on  the  microscale  to  permit  accurate  failure  load  prediction  using 
fracture  mechanics. 

(2)  Limit  analysis  methods  applied  to  microscale  representative  ele¬ 
ments  of  microflawed  materials  permits  the  development  of  limit 
surfaces,  in  general,  three-dimensional  stress  states.  These 
limit  surfaces  are  functions  of  the  flaw  size,  shape  and  spacing 
as  well  as  the  strength  of  the  ductile  matrix  material.  Further, 
these  limit  surfaces  agree  well  with  experimental  data  from  smooth 
specimen  tests. 

(3)  Microscale  limit  surfaces  for  gray  cast  iron,  when  extended  to  the 
macroscale  and  used  in  limit  analysis  of  notched  parts,  provides  an 
upper  bound  to  the  actual  failure  loads.  It  has  been  shown  that  gray 
cast  iron,  especially  in  bending  situations,  is  not  ductile  enough  to 
permit  the  strict  application  of  limit  analysis  tools. 
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5.2  Recommendations  for  Future  Study 

(1)  In  the  case  of  the  RVE  of  the  hole  model,  better  velocity  fields  for  the 
upper  bound  are  needed.  The  object  of  a  study  in  this  area  should  be 
to  cut  off  the  upper  bound  von  Mises  cylinder  in  triaxial  tension  and 
compression  in  order  to  get  the  lower  and  upper  bounds  closer  in  these 
octants. 

(2)  For  real  materials  with  a  distribution  of  flaw  sizes,  an  experimental 
study  is  necessary  to  determine  which  measure  of  flaw  size  distribu¬ 
tion  will  provide  the  best  correlation  with  the  proposed  theory  and 
real  materials. 

(3)  For  the  case  of  the  slit  model,  better  bounds  are  needed.  The  object 
of  study  here  would  be  to  obtain  closer  bounds  in  combined  stress  situ¬ 
ations  involving  tension  and  compression  stresses. 

(4)  An  upper  bound  analysis  is  necessary  for  the  slit  model  when  the  slits 
form  a  spherical  surface  in  the  RVE. 

(5)  Limit  analysis  techniques  should  be  extended  to  cases  involving  adja¬ 
cent  RVE's  with  different  properties.  The  object  here  is  to  model 
materials  which  are  nonhomogeneous.  This  would  enable  the  theory 

to  be  extended  to  quickly  cooled  cast  irons  which  have  dendritic  structure. 

(6)  The  modifications  to  the  limit  analysis  of  bending  situations  should  be 
investigated  more  thoroughly.  A  study  should  be  made  of  the  artificial 
extension  of  the  crack  in  bending  cases  to  determine  if  the  amount  of 
the  extension  is  a  material  dependent  property  and  if  it  is  affected  by 


specimen  geometry. 
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TABLE  1 

FRACTURE  TOUGHNESS  DATA  FOR  DUCTILE  AND 
GRAY  CAST  IRON  FROM  BRANDT  (4) 


Specimen 

Tem^. 

°i' 

K 

max 

(ksiN/Im ) 

V 

(ksiN/in. ) 

Ductile  Iron 

1 

-  50 

--- 

14. 5** 

(ASTM  Grade 

2 

-320 

42.4 

27.0 

60-40-18  from 
Specimen  A536)* 

3 

-320 

42.5 

16.2 

Gray  Iron 

I 

-  50 

37.62 

10. 0** 

2 

-320 

23.9 

23.9 

(ASTM  Class  40)* 

3 

-320 

24.  8 

24.4 

*Classifi cation  is  approximate 

**Based  on  5%  offset  on  P-v  record 
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TENSILE  TEST  DATA  -  PEARLITIC  CAST  IRON 
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TABLE  5 

THREE  POINT  BEND  TEST  RESULTS 
S  =  12,89  FC  Fatigue  Precracked 


B  ■-  1.00 

T  Elevated  Temp. 

Specimen 

No. 

Crack  Length 
(a) 

(inch) 

Fracture 
Load,  kips 

(P) 

Depth 

(w) 

(inch) 

K 

max 

(ksiN/in. ) 

Comments 

3PB-1 

0.60 

18.8 

3.22 

44.0 

3PB-2 

0.60 

18.7 

3.22 

43.  7 

3PB-3 

1.29 

8.5 

3.22 

37.  8 

3PB-4 

1,29 

8.3 

3.22 

36.9 

3PB-FC-1 

0.  65 

15.4 

3.22 

40.0 

3PB-FC-2 

0.  65 

16. 6 

3.22 

43.  1 

3PB-FC-3 

0.  65 

16.9 

3.22 

44.  1 

3PB-FC-4 

0.65 

16.3 

3.22 

42.4 

3PB-T-1 

0.50 

19.4 

3.20 

44,0 

200°  F 

3PB-T-2 

0,50 

19.1 

3.20 

43.3 

200°  F 

3PB-T-3 

0,50 

18.9 

3,20 

42.9 

300°  F 

3PB-T-4 

0.  50 

17.1 

3.20 

38.  8 

300°  F 

3PB-T-5 

0.50 

17.4 

3.20 

39.5 

400°  F 

3PB-T-6 

0.50 

17.8 

3.20 

40,4 

400°  F 
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TABLE  6 

COMPACT  SPECIMEN  TEST  DATA 
W  =  2.0  B  -  1.0 


Specimen 

No. 

Crack  Length 

(a) 

(inch) 

Fracture 

Load,  kips 

(P) 

K 

max 

(ksi\n. 

CT-2 

0.90 

6.0 

35.3 

CT-3 

0.90 

5.15 

30.3 

CT-4 

0.90 

4.75 

28.0 

CT-5 

0.90 

6.36 

37.4 

CT-6 

0.90 

6. 18 

36.4 

CT-7 

0.90 

6.02 

35.4 

CT-8 

0.90 

4.30 

25.3 

TABLE  7 


CIRCUMFERENTIAL  NOTCH  TENSION  DATA 


Specimen 

No. 

Nominal 

Diameter 

D 

(in.) 

Diameter 
at  Notch 
d 

(in.) 

Fracture 
Load,  kips 

(P) 

K 

max 

(ksi  \fln. ) 

CNT-1 

0.50 

0.354 

4.00 

13.6 

CNT-2 

0.50 

0.354 

4.70 

16.0 

CNT-3 

0.50 

0.354 

4.  74 

16.1 

CNT-4 

0.50 

0.354 

5.  60 

19.1 

CNT-7 

2.25 

1.59 

56.0 

19.9 

CNT-8 

2.25 

1.59 

64.0 

22.7 

CNT-9* 

1.00 

0,707 

17.0 

20.4 

CNT-10* 

1.00 

0.707 

18. 1 

21.7 

C  NT-11* 

1.00 

0.707 

18.9 

22.7 

CNT-12* 

0.500 

0.354 

4.96 

16.8 

CNT-13* 

0.500 

0,354 

5.01 

17.0 

CNT-14* 

0.500 

0.354 

5.  84 

19.  3 

^Recalculated  data  from  Drexel  lnstiture  of  Technology 


-  ■ncMMMe  m&mM 


60 


TABLE  8 

DOUBLE  EDGE  CRACK  PLATE  DATA 


B  =  1.00 


Specimen 

No. 

Crack  Length 

(a) 

(inch) 

Fracture 
Load,  kips 

(P) 

Depth 

(w) 

(inch) 

K 

max 

(ksiN/irT. ) 

DEC -2 

0.  75 

36.0 

2.688 

24.6 

DEC-3 

0.  75 

47.0 

2.62 

33.1 

DEC -4 

0.75 

43.9 

2.62 

30.9 

DEC-5 

0.95 

48.5 

3.01 

35.4 

,  i*k& 
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TABLE  9 

PEARLITIC  GRAY  CAST  IRON  DATA  SUMMARY 


Specimen 

No. 

K 

max 

ksi*/in 

y 

ksWin. 

3PB-1 

44.0 

3PB-2 

43.7 

3PB-3 

37.8 

3PB-4 

36.9 

3PB-FC-1 

40.0 

3PB-FC-2 

43. 1 

3PB-FC-3 

44.1 

3PB-FC-4 

42.4 

3PB-T-1 

44.0 

3PB-T-2 

43.3 

3PB-T-3 

42.9 

3PB-T-4 

38.8 

3PB-T-5 

39.5 

3PB-T-6 

40.4 

CT-2 

35.3 

CT-3 

30.3 

CT-4 

28.0 

CT-5 

37.4 

CT-6 

36.4 

CT-7 

35.4 

CT-8 

25.3 

P-P6-1 

21.0 

P-P6-2 

18.  8 

P-P6-3 

32.0 

P-CNT-1 

13.6 

P-CNT-2 

16.0 

P-CNT-3 

16. 1 

P-CNT-4 

19.1 

P-CNT-7 

19.9 

P-CNT-8 

22.7 

P-CNT--9* 

20.4 

P-CNT-10* 

21.7 

P-CNT-11* 

22.7 

P-CNT-12* 

16.  8 

P-CNT-13* 

17.0 

P-CNT-14* 

19.  8 

P-DEC-2 

24.6 

P-DEC-3 

33.1 

P  DEC-4 

30.9 

P-DEC-5 

35.4 

^Recalculated  data  from  Drexel  Institute  of  Technology 
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Fig.  7  Apparent  Microscale  Plastic  Action  at  Graphite  Flake  Tips 
in  Gray  Cast  Iron 

(a)  56%  of  Fracture  Load 

(b)  69%  of  Fracture  Load 

(c)  81%  of  Fracture  Load 
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Fig.  10  Assumed  Stress  Field  for  Lower  Bound  of  Hole  Model 
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(A)  LOWER  BOUND  FOR  HOLE  MODEL 


(B)  PLANE  STRESS  LOWER  AND  UPPER 
BOUND  SURFACE  FOR  HOLE  MODEL 


Fig.  11  Limit  Surfaces  for  RYE  of  Hole  Model 


(A)  DISTRIBUTION  OF 

UNIFORM  SIZE  SPHERES 


© 

© 

0 

© 

d 

p  © 

CROSS  SECTIONAL  VIEW 
OF  A  THREE  DIMENSIONAL 
DISTRIBUTION  OF  UNIFORM 
SIZE  SPHERES 


Fig.  12  Features  oi'  a  Distribution  of  Constant  Radius  Spheres 
in  a  Material 


START  OF  SOLIDIFICATION 
AUSTENITE  DENDRITES  AND 
GRAPHITE  NUCLEI  IN  LIQUID 


AUSTENITE  DENDRITES  AND 
GROWING  EUTECTIC  CELLS 
OF  AUSTENITE  AND  GRAPHITE 


GRAPHITE  FLAKES  IN  MATRIX 
OF  AUSTENITE 


Fig.  13  Schematic  of  Steps  in  the  Solidification  of  Gray  Cast  Iron 


(A)  STRESS  FIELD  IN  THREE  DIMENSIONS 
FOR  LOWER  BOUND  OF  SLIT  MODEL 
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NORMAL1  DISPLACEMENT 
(A)  VELOCITY  FIELD  FOR  ROW  OF  RVE’s 


2  - 


(B)  VELOCITY  FIELD  AT  ONE  RVE  LIGAMENT 


(0  LENGTH  FOR  USE  IN  DISSIPATION 
(MULTIPLY  BY  4) 

Fig.  18  Upper  Bound  Displacement  Field  for  Slit  Model 
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Fig.  21  Plane  Stress  Failure  Surfaces  of  Currently  Available 
Failure  Theories  for  Microflawed  Materials 
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xxx  PROPOSED  a  FISHER 

(PLANE  STRESS) 


Fig.  23  Experimental  Fracture  Data  from  Gray  Cast 
Iron  Study  of  Coffin 
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+  cr2 


Fig.  25  Experimental  Fracture  Data  on  Gray  Iron  from 
Clough  and  Shank 


AXIAL  STRESS  (KSI) 


Fig.  26  Experimental  Fracture  Data  on  Nodular  Iron 
from  Cornet  and  Gras  si,  and  Clough  and  Shank 
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BOUND  OF  DEC  SPECIMEN 
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Fig.  28  Assumed  Stress  and  Velocity  Fields  for 

Limit  Analysis  of  Double  Edge  Crack  Plate 


Fig.  29  Notch  Area  versus  Failure  Load  for  Circumferential  Notch  Tensile  Bars 
Showing  LEFM  and  Proposed  Theory  Predictions 
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K°irt  b~  aT 


(A)  ASSUMED  STRESS  FIELD  FOR  LOWER 
BOUND  OF  3PB  SPECIMEN 


(B)  ASSUMED  VELOCITY  FIELD  FOR  UPPER 
BOUND  OF  3  PB  SPECIMEN 

Fig.  30  Stress  Field  and  Velocity  Field  Assumed  for  Limit 
Analysis  of  a  Notched  Three  Point  Bend  Beam 


CURRENT  STUDY 


Crack  Length  versus  Failure  Load  for  Three  Point  Bend  Tests 
Showing  LEFM  and  Proposed  Theory  Predictions 
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APPENDIX  A 

FUNDAMENTALS  OF  LINEAR  ELASTIC  FRACTURE  MECHANICS 


The  fundamental  concept  in  linear  elastic  fracture  mechanics  is  to  focus  atten¬ 
tion  on  the  material  in  the  vicinity  of  a  crack  tip  because  the  material  behavior  here 
determines  the  material's  resistance  to  fracture.  It  can  be  shown  that  the  elastic 
stress  distribution  at  the  tip  of  a  crack  can  lie  written  as  (42) 


a  - 
x 


a  = 

y 


K  1 

1/2 

2  7r  r  ’  j 

K  1 

1/2 

2  7r  r '  j 

’(a  + 

v  X 

o  ) 
y 

f(0) 

f(6) 


(Al) 

(A2) 

(A3) 


for  the  plane  strain  case.  For  plane  stress, 


(A4) 


In  the  above  equation,  r’  is  the  radial  distance  from  the  crack  tip  to  the  point  of 
interest,  0  is  the  angle  from  the  crack  plane  as  shown  in  Fig.  Al.  The  directions 
x,  y,  and  z  are  also  shown  in  Fig.  Al.  v  is  Poisson's  ratio,  and  K  is  called  the 
stress  intensity  factor.  Tne  K  and  f(0)  values  depe  d  on  the  direction  of  relative 
motion  of  the  crack  surfaces.  If  the  crack  surfaces  open  normal  to  each  other  as 
shown  in  Fig.  A2,  mode  1  crack  displacement  has  occurred.  In  this  case,  K  Is 
given  a  subscript  1.  Kj  is  called  the  mode  I  stress  intensity  factor.  Also  for  mode  I, 
f(0)  becomes 


f(0) 


0 

cos  2 


1  ± 


Restricting  attention  to  the  plane  crack,  0=0°  gives 


(A5) 


f(6)  =  1 


(A  b) 
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Equations  A1  and  A2  now  become 


and 


cr 

x 


K 


I 


(2  7T  r') 


172 


cr 

y 


_Ki 

(2  7rr')1/2 


(A7) 


(A8) 


for  mode  I  crack  surface  displacement. 

From  dimensional  considerations  of  Eqs.  A1  and  A2  or  A7  and  A 8,  the  units 

of  the  stress  intensity  factor  must  be  pounds  per  square  inch  multiplied  by  square 

root  of  inches  assumii  g  that  a  and  a  are  in  pounds  per  square  inch  and  r'  is  in 

x  y 

inches.  In  a  given  piece  of  material  containing  a  crack,  the  most  readily  measured 
parameters  are  the  gross  stress,  a  ,  that  is,  stress  determined  as  if  the  crack 

o 

were  not  present,  and  the  crack  length,  a.  Thus,  the  easiest  representation  for 
the  stress  intensity  factor  is 

K  =  a  a  \Ta  (A9) 

g 


where  a  is  a  geometry  factor  depending  on  mode  of  crack  opening  and  the  loading 
geometry  for  the  specimen  or  part.  It  can  also  be  shown  from  elasticity  ;heory  that 
the  form  of  the  stress  intensity  factor  in  the  general  case  is  as  written  in  Eq.  A9. 

If  a  fixed  geometry  and  crack  size  is  chosen,  a  and  a  become  fixed  in  Eq.  A9. 
The  maximum  gross  stress  that  this  material  can  carry  is  called  the  critical  stress, 
a£.  The  corresponding  stress  intensity  is  called  the  critical  stress  intensity,  Kc< 

Kc  is  a  function  of  material  thickness  as  indicated  in  Fig.  A3.  Above  some  thick¬ 
ness,  K£  is  relatively  constant.  This  constant  is  called  the  plane  strain  fracture 
toughness,  .  Fracture  toughness  is  a  measure  of  the  material's  resistance  to 
fracture  or  rapid  crack  propagation. 
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It  is  generally  accepted  that  Kj  is  a  material  property  and  is  independent  of 
the  geometry  of  loading  or  the  geometry  of  the  crack.  Therefore,  with  the  appro¬ 
priate  a  and  a  value  in  Eq.  A9  and  using  Kj  ,  the  failure  stress  for  any  loading 
geometry  can  be  obtained.  Again,  this  can  be  done  only  if  the  part  being  evaluated  is 
of  sufficient  thickness. 

Strictly  speaking,  linear  elastic  fracture  mechanics  is  not  applicable  in  cases 
where  plastic  zones  appear  at  the  tip  of  a  crack  prior  to  rapid  crack  advance.  How¬ 
ever,  it  has  been  found  that  small  amounts  of  plastic  action  can  be  accounted  for  by 
artificially  extending  the  actual  crack  by  some  amount.  This  amount  is  called  the 
plastic  zone  correction,  r  ,  where,  for  plane  strain 


r 


y 


(A  10) 


In  this  equation  is  the  material's  uniaxial  yield  stress.  The  ry  value  is  utilized 
in  K  calculations  by  replacing  the  crack  length  a  by  a  +  r^.  This  leads  to  an  equa¬ 
tion  involving  K  on  the  left  side  and  K  on  the  right  side  as  indicated  below  with  Kjc 
substituted  for  K. 


a  a 


1 

K 

2 

a  +  w 

a 

- 

ys 

1/2 


(All) 


Utilization  of  this  equation  requires  that  one  guess  at  a  Kjc  value  to  use  on  the  right 
side,  calculate  the  KIc  on  the  left  side,  then  readjust  the  Kjc  on  the  right  until  both 
Kjc  values  arc  equal. 


z 


Fig.  A1  Coordinate  System  at  the  Tip  of  a  Crack 
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APPENDIX  B 

FRACTURE  TOUGHNESS  TESTING 


The  objective  of  fracture  toughness  testing  is  to  determine  the  plane  strain 
fracture  toughness,  K^.  This  objective  is  attained  by  testing  specimens  of  specific 
geometry  to  failure.  Typical  fracture  toughness  specimens  are  shown  in  Fig.  Bl. 
These  specimens  are  usually  cyclically  loaded  at  low  load  prior  to  testing  to  insure 
a  sharp  crack.  If  a  sharp  crack  is  not  present,  the  fracture  toughness  test  measures 
initiation  as  well  as  crack  propagation.  This  is  an  undesirable  situation  when  resist¬ 
ance  to  fracture  (fracture  toughness)  is  to  be  measured. 

In  fracture  toughness  testing,  it  is  necessary  to  measure  the  load  applied  to 
the  specimen  and  the  opening  of  the  crack  surfaces  as  illustrated  for  a  typical  geom¬ 
etry  in  the  sketch  in  Fig.  B2  (43).  The  plot  of  load  versus  crack  surface  displace¬ 
ment  is  called  the  load -displacement  or  P-v  record.  Typical  P-v  plots  are  illus¬ 
trated  in  Fig.  B3  (43).  The  idea  behind  this  record  is  to  enable  the  investigator  to 
determine  the  load  at  which  the  crack  has  extended  2%  of  its  original  length  and  to 
operationally  distinguish  crack  advance  from  plastic  action  (43). 

In  order  to  determine  the  crack  advance  load,  a  secant  line  is  drawn  on  the 
record  as  shown  in  Fig.  B3.  This  secant  has  a  slope  which  is  5%  less  than  the  slope 
of  the  original  linear  portion  of  the  record.  The  load  at  which  the  secant  line  crosses 


the  P-v  plot  is  called  P^.  If  higher  load  appears  on  the  record  prior  to  the  P^  load, 

this  higher  load  is  called  Pq.  If  a  higher  load  is  not  reached  prior  to  the  P&  load, 

then  Pc  -  P„.  The  maximum  load  appearing  on  the  P-v  record  is  P  .  If  P 

b  °  max  max 

occurs  prior  to  P5  then  Pq  =  Pmax.  Pq  is  the  load  used  to  calculate  a  stress  inten¬ 


sity  labeled  Kq.  Kq  is  called  the  "conditional  K^. "  Whether  or  not  Kq  is  a  'valid 
Kjc"  depends  upon  two  criteria.  First,  Pmax/PQ  must  be  less  than  1.10.  Secondly, 
both  the  thickness  and  crack  length  of  the  specimen  must  be  greater  than  2.  5  (Kq/ct  )' 
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If  these  criteria  are  met,  then  Kq  -  K^  (43).  The  purpose  of  this  operation!  1 
exercise  with  the  P-v  record  is  to  insure  that  the  specimen  tested  was  of  sufficient 
si/.e  to  give  a  plane  strain  fracture  toughness  and  not  some  other  artificially  higher 
toughness.  The  terminology  'valid  K^  s  is  used  to  indicate  that  plane  strain  condi¬ 
tions  have  indeed  been  met. 

The  preceding  discussion  oi  P-v  record  is  according  to  an  ASTM  criteria 

published  in  1972  (43).  The  previous  (1969)  criteria  differs  slightly  from  the  later 

standard.  The  criteria  for  valid  in  1969  was  that  the  horizontal  distance  between 

the  actual  record  and  the  secant  corrected  line  at  80%  of  P,.  must  be  less  than  25% 

of  the  corresponding  distance  at  the  P5  load.  This  procedure  is  illustrated  in  Fig.  B4. 

This  procedure  was  replaced  in  1972  by  the  Pm  /P„  criteria.  In  this  study  the  ear- 

max  u  j 

lier  1969  criteria  will  be  used  throughout. 

In  general,  if  the  requirements  of  the  above  secant  correction  procedure  are 
not  met  and  Kq  is  not  equal  to  K^,  then  the  plastic  action  occurring  at  the  tip  of 
the  crack  involves  a  considerable  volume  and  the  equation  and  assumptions  of  linear 
elastic  fracture  mechanics  are  not  applicable. 


(A)  CIRCUMFERENTIAL 
NOTCH  TENSION 


(B)  DOUBLE  EDGE  CRACK 
PLATE  (DEC) 


(C)  THREE  POINT  BEND  (3PB) 


(D)  COMPACT  SPECIMEN  (CT) 

Fig.  B1  Typical  Fracture  Toughness  Test  Specimens 


Fig.  B3  TVpical  P-v  Records 


